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ABSTRACT 


In this paper we have analysed further some 
of our previously published data relating to 
thermal effects produced on the ignition of a 
range of pyrotechnic compositions, and have 
evaluated the hazards posed to those handling 
and working with such materials by reference 
to the distances for different degrees of burn 


injury. 


Keywords: thermal hazard, radiation, burn 
injury, pyrotechnics 


Introduction 


The thermal characteristics of fireballs from 
a number of fuel sources other than pyrotech- 
nics have been reported in the literature, includ- 
ing liquid propellants,''! motor fuels'”! and pro- 
pane.!! Fireball effects can be described in 
terms of the maximum diameter D (m), the du- 
ration of the thermal effect ¢ (s) and the heat 
flux radiated from the fireball surface g (kW m- 
*). In many cases both D and t are related to the 
fuel mass M (kg) by a power law relationship of 
the approximate form M' 


To investigate the applicability of this gen- 
eralised equation to pyrotechnics, a series of 
experiments was recently done! with a range 
of compositions in which D, ¢ and g were quanti- 
fied at various values of M, up to 25 kg. The 
results from these trials indicated that relation- 
ships of the form D = aM“ and t= bM* applied. 
However, the values of x varied from 1/3, 
sometimes significantly, in both equations, and 
different compositions gave different values for 
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the constants a and 6 (within the range 0.53— 
26.0). 


The predictive equations obtained from the 
study were used in a subsequent paper to exam- 
ine both the fire and explosion hazards of the 
pyrotechnics that had been tested.'*! Explosion 
hazards were evaluated since some of the com- 
positions exploded under conditions of self con- 
finement. The Eisenberg thermal radiation dose 
criterion'®! was used to evaluate potential levels 
of harm (blister thresholds and degrees of burn) 
at different distances. 


In this short paper, some of the implications 
of these results are examined with respect to 
those in the pyrotechnics industry working with 
relatively small quantities of loose (self con- 
fined) composition. 


Discussion 


Table 1 lists the compositions of the pyro- 
technics examined. The potential hazards of 
pyrotechnic fireballs can be defined in terms of 
heat flux, time of exposure (or burning time, if 
less) and distance from the fireball surface since 
distance determines thermal exposure and dura- 
tion influences the dose received. 


The various design requirements of pyro- 
technic mixes are reflected partly in their burn 
times and these show a significant range in du- 
ration from ¢ <1 s to >50 s for 1-25 kg quanti- 
ties of different materials.'*! 


It has been reported in the literature'”! that it 
takes approximately 5 s to sense high levels of 
thermal radiation and start to make an escape, 
and it is also known that a proportion of the 
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Table 1. List of Pyrotechnic Compositions and Their Ingredients. 


Pyrotechnic Substance Ingredient % by Mass 


Gunpowder 3/7 Grist 
Carbon 
Sulfur 
Flare Composition 1 


Flare Composition 2 


Potassium nitrate 


Sodium nitrate 
Calcium oxalate 


Star Composition 1 


Potassium perchlorate 


PVC powder 


Star Composition 2 


Aluminium 
Potassium nitrate 


Priming Composition 1 


Silicon powder 
Gunpowder sulfurless mealed 


Priming Composition 2 


Gunpowder 


Potassium nitrate 


Silicon 


Dextrin binder 


Flash Composition 1 


Potassium perchlorate 


Graphite 


final burn injury can occur during the phase 
when the skin is cooling.'”! 


It may be possible for workers exposed to 
the thermal radiation from fireballs with t >5 s 
to limit their potential thermal dose by making 
a rapid escape. On the other hand, certain pyro- 
technic materials burn for less than 5 s and for 
such events an eventual escape would not be 
likely to limit the received thermal radiation 
dose. 


To compare the potential of different pyro- 
technics to cause burn injuries, with fireballs 
burning for greater than 5 s we have still used 
the 5s reaction time!” in our calculations since, 
although a process operator seated directly in 
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front of an igniting pyrotechnic material is 
likely to respond in less than 5 s, the overall dura- 
tion of exposure can be assumed to be ap- 
proximately 5 s. The extreme conditions arising 
from engulfment within the fireball have been 
assumed to be fatal. 
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Flash Comp. 1 


Star Comp. 1 
Gunpowder 


Priming Comp. 2 
Flare Comp. 2 
Star Comp. 2 


Flare Comp. 1 


Priming Comp. 1 


2.5m 


Distance from Centre of Fireball (m) 


Figure 1. Flare / fire / fireball dimensions (m) for 1 kg quantities of loose pyrotechnic compositions. 


From our previously published results," it is 
possible to compare the sizes of actual and 
pseudo fireballs produced by the ignition of dif- 
ferent compositions with M = 1 kg, Figure 1. 
This diagram indicates that most of the materi- 
als examined produce fireballs that would en- 
gulf a process worker positioned at arms length 
from the point of initiation. Dose levels and cor- 
responding burns injuries for process workers 
positioned beyond the fireball radius can be cal- 
culated using a maximum exposure time of 5s. 
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A dosage of 1200 (kW m7”)*s has been used in 
the literature!” as the mean value to produce 
second degree burns with depths >0.1 mm on 
unprotected skin and this value was employed 
in our work. Clearly, appropriate fire protective 
clothing can offer some mitigation, but neverthe- 
less, this dose has been reported to result in 1% 
lethality to averagely dressed exposees.'” The 
distances to second degree burns for every 
composition with M = | kg, Figure 2, show that 
other process workers in the same room could 
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Star Comp. 1 


Flare Comp. 2 


Priming Comp. 2 
Star Comp. 2 


Gunpowder 


Flare Comp. 1 
Priming Comp. 1 


0 1 


1.7 |21 3.1 


4.3m 


2 3 4 


ee eee eee eee eee 
Distance from Centre of Fireball (m) 


Figure 2. Distances (m) to second degree burns from 1 kg quantities of loose pyrotechnic 


compositions. 


well be within the minimum distances to be 
affected. The information summarised in Figure 
2 should, however, be taken as a first approxi- 
mation since the hazards posed by the slowest 
burning materials may be overestimated because 
the flames can take time to build up to peak ir- 
radiance (i.e., the calculations assume 5 s of ex- 
posure to the peak surface emissive power from 
fireballs). 


Since the blister threshold dosage'” is only 
210(kW m”)*"s, such injuries can be sustained 
at greater distances from the point of initiation 
of a range of substances with M = 1 kg, Figure 
3. Again, other workers in the same process 
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room and at considerable distances from the 
source could receive blister injuries. 


The importance of taking account of exposure 
time in determining dosage is illustrated by the 
different rankings obtained for the potential 
hazards posed by the fireballs from 1 kg quanti- 
ties of materials when using either distance to 
second degree burns (dosage) or fireball dimen- 
sions as the ranking criterion, Table 2. This Ta- 
ble shows that a simplistic ranking based on 
fireball dimensions alone may not accurately 
represent the hazard posed to workers outside 
the fireball diameter by the ignition of certain 
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Star Comp. 1 


Flare Comp. 2 


Priming Comp. 2 
Star Comp. 2 
Gunpowder 


Flare Comp. 1 
Priming Comp. 1 


S30 ‘a6 5.9 8.3 
3.8 
3 4 5 6 67 8 


Distance from Centre of Fireball (m) 


Figure 3. Distances (m) to blister threshold for I kg quantities of loose pyrotechnic compositions. 


pyrotechnic compositions (e.g., gunpowder and 
flare composition 2). 


By assuming that our fireball data for gun- 
powder with M = | to 25 kg can be extrapolated 
to lower values of M, it is possible to use the 
relationship"! D = 3.1 M°?” to estimate fireball 
sizes and hence the hazards posed by the igni- 
tion of small quantities of gunpowder, as indi- 
cated in Figure 4. This diagram shows that a 
fireball produced by burning quantities as low 
as 100 g can engulf a process worker. 


In terms of reducing potential hazards in the 
working environment posed by the accidental 
ignition of pyrotechnics, the first step is to 
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minimise the quantities being handled. This is 
often supplemented by the introduction of engi- 
neering controls to ensure that a suitable physi- 
cal separation exists between the worker and 
the potential heat source. Techniques employed 
include remote handling, the use of robotics, 
and the use of protective screens. Some fire pro- 
tection can also be gained by using quenching 
systems, 8.191 


As a last resort, personal protective equip- 
ment can be employed and a recently published 
guide gives an indication of the current position 
with regard to the selection and use of fire pro- 
tective clothing for explosives workers in the 
UK.!"'"! The ICARUS code!'”! is also able to 
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Table 2. Relative Rankings of the Hazards Posed by Fireballs from 1 kg Quantities of Pyrotech- 


nic Compositions. 


Ranking Based 
on Fireball 
Dimensions 


Fireball 


Pyrotechnic 
Star comp. 1 


Flare comp. 2 
Star comp. 2 
Flare comp. 1 
Priming comp. 1 


predict burn injuries for workers exposed to 
thermal radiation when wearing appropriate 
protective clothing. 


Research sponsored by the Health and Safety 
Executive is currently underway to develop a 
full torso portable manikin to enable the evalua- 
tion of complete garments against the thermal 
threat posed by burning pyrotechnics and pro- 
pellants, and this may provide a useful means of 
ranking the performance of protective clothing, 
thus aiding selection. 


Conclusions 


In this paper the potential thermal hazards 
posed to workers in the pyrotechnics industry 
by the ignition of different types of composi- 
tions have been examined. Relatively small 
quantities of material can produce a significant 
fireball and, because thermal emissive powers 
are relatively large, even short duration expo- 
sure can result in burn injuries at considerable 
distances from the source. 


The main factors affecting hazard are the 
close proximity of the pyrotechnic worker to 
the composition and the quantity present in the 
process room. Means of minimising the role of 
these factors (i.e., reducing quantity, increasing 
distance) are desirable in terms of improving 
safety in the workplace. 
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Distance to Sec- 
ond Degree 


Ranking Based 
on Level of Harm 
(Second Degree 

Burns 
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Errata — Number 5 


Page 3, Figure 3 Incorrect subscript for “h’’. Is should have been “con”. 


Figure 3. Cross section of a rocket motor with a conical combustion cavity. 
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Methods for Calculation of Thrust Coefficients 


Ed Brown 
PO Box 177, Rockvale, CO 81244, USA 


ABSTRACT 


When one wishes to predict or analyze rocket 
motor/engine performance accurately, a method 
for determining the nozzle thrust coefficient at 
each point in time is essential. If the vacuum 
thrust coefficient for the nozzle is known, the 
thrust coefficient for other conditions can easily 
be determined. This paper outlines several 
methods for determining the vacuum thrust coef- 
ficient in the hope that it will encourage this 
type of calculation and lead to new and better 
methods for simulation and analysis in the fu- 
ture. 


Key Words: vacuum thrust coefficient, thrust 
coefficient, expansion ratio, specific heat ratio, 
chamber pressure, thrust, rocket propulsion the- 
ory 


Introduction 


This paper will neither present nor attempt to 
explain any of the theory involved, as refer- 
ences''”! exist that do this more than ade- 
quately and far better than this writer could. It 
will attempt to outline and present methods that 
may be used to simplify basic design and analysis 
calculations. While graphs"! and tables"! exist 
for thrust coefficients etc., determining more ac- 
curate values than can be obtained from a graph 
or table is often desirable. Also, incorporating 
their real time calculation in a simulation or 
analysis program or procedure is often benefi- 
cial, and graphs and tables are not always avail- 
able or convenient to use. One equation'*! used 
to calculate values for the thrust coefficient (Cr) 
is: 
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where capital gamma (I), also known as the 
Vandenkerckhove function,'! is defined by: 


(2) 


cy 
2 [R2(r-) 
ills 
in addition, gamma (y) is the ratio of specific 
heats (or specific heat ratio) for the exhaust 
gases. The pressure at the nozzle exit plane (P.), 
the chamber pressure (P,), and the ambient pres- 
sure (P,) need to be in absolute pressure units. 
Nozzle exit area (A,) and nozzle throat area (A, ) 
should be given using the same units for each. 
Nozzle expansion ratio”! () can be defined as: 


€=—* (3) 


(4) 


Equation 4 can be shortened by defining alpha 
(C) as: 


Oo = 7 (5) 
beta (f) as: 
p=1— (6) 
Hf 


and delta (A) as: 
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A= — (7) 


(8) 


Except for pressure terms, all variables in equa- 
tion 8 are defined in terms of the specific heat 
ratio (y). The specific heat ratio is usually known 
or can be calculated using one of the many PEP 
programs (PROPEP, CETPC, ISP, SP-273 etc.). 
However, if it is not known, a value of 1.25 can 
be used (not a bad assumption for most solid 
propellants used in high power rocketry). 


If atmospheric pressure is reduced to zero in 
equation 1, then the thrust coefficient becomes 
the vacuum thrust coefficient (C;, _): 


This allows equation | to be rewritten in a 
simpler form: 


P, 
C=C, > {=| (10) 


Tools are now in place for further calculations. 


Short BASIC program listings are included 
in the appendices to illustrate each calculation 
method described. The simplest, a method for 
calculating the optimum expansion ratio and 
thrust coefficient for a fixed pressure ratio is 
outlined first. A method for calculating thrust 
coefficients using a fixed expansion ratio is 
then presented. Following this is an approxi- 
mate method for calculation of thrust coeffi- 
cients (again using a fixed expansion ratio). The 
final method is for converting thrust levels to 
chamber pressures. This method can be re- 
versed to predict thrust versus time in simula- 
tion programs (after using Kn values to deter- 
mine pressures). For all methods, it is assumed 
that a specific heat ratio has been provided for 
use. 
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A Method for Calculation of Opti- 
mum Thrust Coefficient Using a 
Fixed Pressure Ratio!pren™* 4 


The easiest thrust coefficient calculation is 
determination of the optimum expansion ratio 
and thrust coefficients for a fixed pressure ratio. 
Ambient pressure at the operation or simulation 
location of the motor must be known or se- 
lected. Motor operating pressure is then se- 
lected. A nozzle’s optimum thrust coefficient is 
obtained when its exit pressure is equal to am- 
bient pressure. The pressure ratio is calculated 
using: 


P,P 
PRe oS ee (11) 
e Le 


Values for equations 5, 6, and 7 are determined 
and the results, to this point, are inserted into 
equations 2 and 8 to produce values for capital 
gamma and the expansion ratio. Final steps are 
the calculation of the vacuum thrust coefficient 
and thrust coefficient using equations 9 and 10. 


A Method for Calculation of 
Thrust Coefficient Using a 
Specific Expansion Ratio “Pren™* 8 


In the real world, the expansion ratio is usu- 
ally a fixed value for initial design and testing, 
with possible optimization occurring during final 
development. One very accurate method of cal- 
culating vacuum thrust coefficients is to select 
an expansion ratio and then calculate a tempo- 
rary expansion ratio using an arbitrarily chosen 
pressure ratio. This result is then compared with 
the selected expansion ratio. If they do not 
match, a new temporary expansion ratio is cal- 
culated using a new pressure ratio (selected to 
ensure convergence of expansion ratio values). 
This process is repeated until the temporary 
expansion ratio matches the selected expansion 
ratio. This pressure ratio is then used to calcu- 
late the vacuum thrust coefficient using equa- 
tion 9. The thrust coefficient for any chamber 
operating pressure can now be calculated using 
equation 10, if ambient pressure is known for 
motor operation location. This method gives 
results for vacuum thrust coefficients and pres- 
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sure ratios that have very good agreement with 
tables."#! 


A Brute Force Method for 
Calculation of Thrust Coefficients 
Using Specified Expansion 
Ratios APrendix Cj 


The writer believes that calculating the pres- 
sure ratio directly from the specific heat ratio 
and expansion ratio may be advantageous. 
While this may not be true, he has been unable 
to accomplish this to date. An attempt to “force” 
a solution to the problem has been made. Using 
the method outlined in Appendix B, vacuum 
thrust coefficients and pressure ratios were cal- 
culated for expansion ratios ranging from one to 
eight (by increments of 0.5) for specific heat 
ratios ranging from 1.1 to 1.4 (by increments of 
0.5). A shareware program, CurveFit'” was 
then used for each ratio of specific heat to relate 
pressure ratio and expansion ratio. An equation 
having the form: 


PR=a+be+ce* (12) 


seemed to fit each group of data. A new prob- 
lem was that a, b, and c had differing values for 
each specific heat ratio. In an attempt to make 
the results more universal, another round of 


curve fitting was attempted to relate each coef- 
ficient to specific heat ratios. This resulted in 
the following results: 


_ 38.7019 poe! See (13) 
7 % 
= _ (14) 
0.3283(y — 1.5242)” + 0.1423 
c = (52.4092)y'** (0.0013) " (15) 


The equations for coefficients a and b were the 
first listed by the program while the equation 
for c was the third listed (but appeared to fit the 
data almost as well and was much easier to cal- 
culate). 


Next, a short program was written using 
these equations to see if results were compara- 
ble to table values or those produced using the 
method in Appendix B. The program produced 
very good agreement for vacuum thrust coeffi- 
cients, but only good agreement for pressure 
ratios. Since these pressure ratios are only being 
used for calculating vacuum thrust coefficients, 
this method may be an acceptable substitute for 
other calculation methods. Table 1 compares re- 
sults from programs in Appendix B and Appen- 
dix C using a specific heat ratio of 1.25. 


Table 1. Comparison of Results Using Appendix B and Appendix C Methods. 


Bs Pa B 


C Appendix B a ee c Appendix C 
(Cr. (PR) Chin. (PR) 


1.248592 
1.391232 
1.463144 


1. "680943 
1.688710 
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1.680943 
1.688710 
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Determining Operating Pressures 
from Thrust Levels Using Vacuum 
Thrust Coefficients /* 4Pren* PI 


Thrust levels used to calculate chamber 
pressures may come from thrust curve data, 
manufacturer supplied information or any other 
reasonable source. Nozzle throat and exit areas 
are needed (easily calculated if we measure or 
are given throat and exit diameters). A simple 
equation”! for thrust (F) is: 


F=C,P.A, (16) 


Rearranging to solve for chamber pressure 
gives: 
Y ioe Z (17) 
“ CrA, 
Multiplying both sides of equation 10 by P, 
gives: 


C,P, =C, P. -8P, (18) 


Rearranging this to solve for chamber pressure 
(P..): 
p= Sect (19) 
Cr. 
Reexamining equation 17 and multiplying 
both sides by Cr produces: 
F 


CP, = (20) 


t 


which can easily be solved using throat areas 
and thrust levels. Inserting this and earlier re- 
sults into equation 19 produces a value for 
chamber pressure. While many other factors 
can (and perhaps should) be considered,!'?! 
chamber pressure values calculated using this 
method give surprisingly good agreement with 
measured pressures. In this writer’s opinion, 
thrust levels are much easier to obtain than 
pressure measurements. 


Summary 


Several methods of calculating vacuum thrust 
coefficients and a method for calculating cham- 
ber pressure have been presented. Once the vac- 
uum thrust coefficient has been calculated, it is 


an easy matter to calculate the thrust coefficient 
for any chamber pressure using equation 10. A 
value for the thrust coefficient is essential for 
simulating or analyzing motor performance. A 
word of warning, if the nozzle is severely over- 
expanded, separation of flow will occur in the 
nozzle with an effective reduction in expansion 
ratio. Although most of the discussion has im- 
plied usage with solid motors, these methods 
apply to hybrid motors or liquid rocket engines 
as well. The programs in the appendices include 
a calculation for approximate characteristic ex- 
haust velocity (c*) if the chamber temperature 
and molecular weight of the exhaust gases are 
known. If not, enter any positive number when 
prompted for them. They do not affect the calcu- 
lation of the other values. The programs are writ- 
ten anticipating the use of English units. It is 
hoped that encouragement of development of 
more and better methods for simulating and 
analyzing motor performance has been accom- 
plished. 
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Appendix A 


10 CLS ' PROGRAM NAME IS APPENDXA.BAS 

20 PRINT "PLEASE ENTER RATIO OF SPECIFIC HEATS." 

30 INPUT "ENTER NUMBER BETWEEN 1.05 AND 1.71"; GAMMA 

40 IF GAMMA >= 1.05 AND GAMMA <= 1.71 THEN GOTO 50 ELSE BEEP: GOTO 20 


50 ' Ae CALCULATE VALUE FOR CAPITAL GAMMA Ae 
60 CAPITAL.GAMMA = SQR(GAMMA) * (2 / (GAMMA + 1)) *~ ((GAMMA + 1) / (2 * 
(GAMMA - 1))) 


70 ALPHA SQR((GAMMA - 1) / (2 * GAMMA)): BETA = (GAMMA - 1) / GAMMA: 
DELTA = 1 / GAMMA 


80 INPUT "ENTER ATMOSPHERIC PRESSURE (PSIA)"; ATMOS.PRESS 

90 INPUT "ENTER CHAMBER PRESSURE (PSIA)"; CHAMBER. PRESS 

100 PRESS.RATIO CHAMBER.PRESS / ATMOS.PRESS 

110 PRESS.RATIO = INT(PRESS.RATIO * 1000 + .5) / 1000 

120 NUMERATOR = CAPITAL.GAMMA * ALPHA * (PRESS.RATIO) * DELTA 

130 DENOMINATOR = SQR(1 - (1 / PRESS.RATIO) * BETA) 

140 =! EI CALCULATE EXPANSION RATIO Oe 
150 EXPANSION.RATIO NUMERATOR / DENOMINATOR 

160 EXPANSION.RATIO = INT(EXPANSION.RATIO * 1000 + .5) / 1000 

170 PRINT "OPTIMUM EXPANSION RATIO IS"; EXPANSION.RATIO 


180 INPUT "ENTER ADIABATIC FLAME TEMPERATURE IN DEGREES RANKINE"; 
FLAME . TEMP 


190 INPUT "ENTER MOLECULAR WEIGHT OF COMBUSTION PRODUCTS"; MOL.WEIGHT: 
PRINT 


200 '! Ae CALCULATE VACUUM THRUST COEFFICIENT aCe 
210 DENOMINATOR = SQR(1 - (1 / PRESS.RATIO) * ((GAMMA - 1) / GAMMA) ) 


220 VACTHRUST.COEFF = CAPITAL.GAMMA * SQR(2 * GAMMA / (GAMMA - 1)) * DE- 
NOMINATOR + EXPANSION.RATIO * 
(1 / PRESS.RATIO) 


230 PRINT "VACUUM THRUST COEFFICIENT IS"; VACTHRUST.COEFF 

240 =! Ee CALCULATE CHARACTERISTIC EXHAUST VELOCITY aCe 
250 C.STAR = SQR(49800! * FLAME.TEMP / MOL.WEIGHT) / CAPITAL.GAMMA 

260 PRINT "CHARACTERISTIC EXHAUST VELOCITY IS"; C.STAR; "FEET PER SECOND." 
270 PRINT "PRESSURE RATIO IS"; PRESS.RATIO: PRINT 

280 THRUST.COEFF = VACTHRUST.COEFF - EXPANSION.RATIO * (1 / PRESS.RATIO) 
290 PRINT "THRUST COEFFICIENT AT"; CHAMBER.PRESS; "PSIA IS"; THRUST.COEFF 
300 INPUT "DO YOU WISH TO CONTINUE (Y/N)"; ANSWERS 

310 IF ANSWERS = "y" OR ANSWERS = "Y" THEN GOTO 10 ELSE END 
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Appendix B 


CLS ' PROGRAM NAME IS APPENDXB.BAS 
PRINT "PLEASE ENTER RATIO OF SPECIFIC HEATS." 

INPUT "ENTER NUMBER BETWEEN 1.05 AND 1.71"; GAMMA 

IF GAMMA >= 1.05 AND GAMMA <= 1.71 THEN GOTO 50 ELSE BEEP: GOTO 20 
wk CALCULATE VALUE FOR CAPITAL GAMMA wk 
CAPITAL.GAMMA = SQR(GAMMA) * (2 / (GAMMA + 1)) * ((GAMMA +1) / (2 * 
(GAMMA - 1))) ! 


ALPHA = SQR((GAMMA - 1) / (2 * GAMMA)): BETA = (GAMMA - 1) / GAMMA: 
DELTA = 1 / GAMMA 
' eee INITIALIZE CHAMBER PRESSURE/EXIT PRESSURE Ae 


LET PRESS.RATIO = 1.72 
PRINT "ENTER SELECTED EXPANSION RATIO." 
INPUT "VALUE MUST BE BETWEEN 1.00 AND 50.00"; EXPANSION.RATIO: PRINT 
: PRINT 
EXPANSION.RATIO = INT(EXPANSION.RATIO * 1000 + .5) / 1000 
IF EXPANSION.RATIO >= 1 AND EXPANSION.RATIO <= 50 THEN GOTO 140 ELSE 
BEEP: GOTO 100 
LET TEMP.EXP.RATIO = 0: COUNT = 0 
; Ae LOOP TO MATCH SELECTED EXPANSION RATIO Ae 
WHILE TEMP.EXP.RATIO <> EXPANSION.RATIO 
IF TEMP.EXP.RATIO < EXPANSION.RATIO THEN PRESS.RATIO = PRESS.RATIO 
+ (EXPANSION.RATIO - TEMP.EXP.RATIO) * 4! + .001 ELSE GOTO 180 
PRESS.RATIO = PRESS.RATIO - (TEMP.EXP.RATIO - EXPANSION.RATIO) * 
4! - .002 
PRESS.RATIO = INT(PRESS.RATIO * 1000 + .5) / 1000 
NUMERATOR = CAPITAL.GAMMA * (PRESS.RATIO) * DELTA * ALPHA 
DENOMINATOR = SQR(1 - (1 / PRESS.RATIO) “* BETA) 
TEMP.EXP.RATIO = NUMERATOR / DENOMINATOR 
TEMP.EXP.RATIO = INT(TEMP.EXP.RATIO * 1000 + .5) / 1000 
COUNT = COUNT + 1 
WEND 
: basta the) CONTINUE WITH FURTHER CALCULATIONS Dae 
INPUT "ENTER LOCAL ATMOSPHERIC PRESSURE (14.7 IF UNSURE)"; ATMOS.PRESS 
INPUT "ENTER ADIABATIC FLAME TEMPERATURE IN DEGREES RANKINE"; 
FLAME . TEMP 
INPUT "ENTER MOLECULAR WEIGHT OF COMBUSTION PRODUCTS"; MOL.WEIGHT: 
PRINT 
’ ste) CALCULATE VACUUM THRUST COEFFICIENT Pe 
VACTHRUST.COEFF = CAPITAL.GAMMA * SQR(2 / BETA) * DENOMINATOR + EXPAN- 
SION.RATIO * (1 / PRESS.RATIO) 
PRINT "VACUUM THRUST COEFFICIENT IS"; VACTHRUST. COEFF 
, aCe CALCULATE CHARACTERISTIC EXHAUST VELOCITY Ae 
C.STAR = SQR(49800! * FLAME.TEMP / MOL.WEIGHT) / CAPITAL.GAMMA 
PRINT "CHARACTERISTIC EXHAUST VELOCITY IS"; C.STAR; "FEET PER SECOND." 
PRINT "PRESSURE RATIO IS"; PRESS.RATIO: PRINT 
INPUT "WHAT IS THE DESIRED CHAMBER PRESSURE"; CHAMBER. PRESS 
THRUST.COEFF = VACTHRUST.COEFF - EXPANSION.RATIO * (1 / PRESS.RATIO) 
PRINT "THE THRUST COEFFICIENT AT"; CHAMBER.PRESS; "PSIA IS"; 
THRUST. COEFF 
INPUT "DO YOU WISH TO CONTINUE (Y/N)"; ANSWERS 
IF ANSWERS = "Y" OR ANSWERS = "y" THEN GOTO 10 ELSE END 
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Appendix C 


CLS ' PROGRAM NAME IS APPENDXC.BAS 

PRINT "PLEASE ENTER RATIO OF SPECIFIC HEATS." 

INPUT "ENTER NUMBER BETWEEN 1.05 AND 1.71"; GAMMA 

IF GAMMA >= 1.05 AND GAMMA <= 1.71 THEN GOTO 50 ELSE BEEP: GOTO 20 

: Ka CALCULATE VALUE FOR CAPITAL GAMMA ka 
CAPITAL.GAMMA = SQR(GAMMA) * (2 / (GAMMA + 1)) * ((GAMMA +1) / (2 * 
(GAMMA - 1))) ! 

BETA = (GAMMA - 1) / GAMMA 

' kK INITIALIZE CHAMBER PRESSURE/EXIT PRESSURE eK 
PRINT "ENTER SELECTED EXPANSION RATIO." 

INPUT "VALUE MUST BE BETWEEN 1.00 AND 50.00"; EXPANSION.RATIO: PRINT 
: PRINT 

IF EXPANSION.RATIO >= 1 AND EXPANSION.RATIO <= 50 THEN GOTO 120 ELSE 
BEEP: GOTO 90 

' kK CALCULATE COEFFICIENTS A, B, AND C kK 
A 38.7019 / GAMMA - 15.5227 / GAMMA * 2 - 25.6441 
B= 1/ (.3283 * (GAMMA - 1.5242) * 2 + .1423) 
C 

1 


= 52.4092 * GAMMA * 15.0965 * .0013 “~ GAMMA 
RICK CALCULATE PRESSURE RATIO Ex 
PRESS.RATIO = A + B * EXPANSION.RATIO + C * EXPANSION.RATIO “* 2: 
PRESS.RATIO = INT(PRESS.RATIO * 1000 + .5) / 1000 


DENOMINATOR = SQR(1 - (1 / PRESS.RATIO) * BETA) 
‘ EER CONTINUE WITH FURTHER CALCULATIONS BK, 
INPUT "ENTER LOCAL ATMOSPHERIC PRESSURE (14.7 IF UNSURE)"; ATMOS.PRESS 


INPUT "ENTER ADIABATIC FLAME TEMPERATURE IN DEGREES RANKINE"; 
FLAME . TEMP 
INPUT "ENTER MOLECULAR WEIGHT OF COMBUSTION PRODUCTS"; MOL.WEIGHT: 
PRINT 
: AX CALCULATE VACUUM THRUST COEFFICIENT ia 
VACTHRUST.COEFF = CAPITAL.GAMMA * SQR(2 / BETA) * DENOMINATOR + EXPAN- 
SION.RATIO * (1 / PRESS.RATIO) 

PRINT "VACUUM THRUST COEFFICIENT IS"; VACTHRUST. COEFF 

: AS CALCULATE CHARACTERISTIC EXHAUST VELOCITY ROE, 
C.STAR = SQR(49800! * FLAME.TEMP / MOL.WEIGHT) / CAPITAL.GAMMA 

PRINT "CHARACTERISTIC EXHAUST VELOCITY IS"; C.STAR; "FEET PER SECOND." 
PRINT "PRESSURE RATIO IS"; PRESS.RATIO: PRINT 

INPUT "ENTER DESIRED CHAMBER PRESSURE"; CHAMBER. PRESS 

THRUST.COEFF = VACTHRUST.COEFF - EXPANSION.RATIO * (AT- 

MOS . PRESS/ CHAMBER . PRESS) 

PRINT "THE THRUST COEFFICIENT AT"; CHAMBER.PRESS; "PSIA IS"; 
THRUST . COEFF 

INPUT "DO YOU WISH TO CONTINUE? (Y/N)"; ANSWERS 

IF ANSWERS = "y" OR ANSWERS = "Y" GOTO 10 ELSE END 
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Appendix D 


10 CLS ' PROGRAM NAME IS APPENDXD.BAS 

20 PRINT "PLEASE ENTER RATIO OF SPECIFIC HEATS." 

30 INPUT "ENTER NUMBER BETWEEN 1.05 AND 1.71"; GAMMA 

40 IF GAMMA >= 1.05 AND GAMMA <= 1.71 THEN GOTO 50 ELSE BEEP: GOTO 20 

50 ' eK CALCULATE VALUE FOR CAPITAL GAMMA kK 

60 CAPITAL.GAMMA = SOR(GAMMA) * (2 / (GAMMA + 1)) * ((GAMMA +1) / (2 * 

(GAMMA - 1))) '! 

70 ALPHA = SOR((GAMMA - 1) / (2 * GAMMA)): BETA = (GAMMA - 1) / GAMMA: 

DELTA = 1 / GAMMA 

80 ' kKK INITIALIZE CHAMBER PRESSURE/EXIT PRESSURE aK 

90 LET PRESS.RATIO = 1.72 

100 PRINT "ENTER SELECTED EXPANSION RATIO." 

110 INPUT "VALUE MUST BE BETWEEN 1.00 AND 50.00"; EXPANSION.RATIO: PRINT 

: PRINT 

120 EXPANSION.RATIO = INT(EXPANSION.RATIO * 1000 + .5) / 1000 

130 IF EXPANSION.RATIO >= 1 AND EXPANSION.RATIO <= 50 THEN GOTO 140 ELSE 

BEEP: GOTO 100 

140 LET TEMP.EXP.RATIO = 0: COUNT = 0 

150 '! KkK LOOP TO MATCH SELECTED EXPANSION RATIO kK 

160 WHILE TEMP.EXP.RATIO <> EXPANSION. RATIO 

170 IF TEMP.EXP.RATIO < EXPANSION.RATIO THEN PRESS.RATIO = PRESS.RATIO 
+ (EXPANSION.RATIO - TEMP.EXP.RATIO) * 4! + .001 ELSE GOTO 180 

180 PRESS.RATIO = PRESS.RATIO - (TEMP.EXP.RATIO - EXPANSION.RATIO) * 

4! - .002 

190 PRESS.RATIO = INT(PRESS.RATIO * 1000 + .5) / 1000 

200 NUMERATOR = CAPITAL.GAMMA * (PRESS.RATIO) * DELTA * ALPHA 

220 DENOMINATOR = SQR(1 - (1 / PRESS.RATIO) * BETA) 

230 TEMP.EXP.RATIO = NUMERATOR / DENOMINATOR 

240 TEMP.EXP.RATIO = INT(TEMP.EXP.RATIO * 1000 + .5) / 1000 

250 COUNT = COUNT + 1 

260 WEND 

270 =! KKK CONTINUE WITH FURTHER CALCULATIONS KKK 

280 INPUT "ENTER LOCAL ATMOSPHERIC PRESSURE (14.7 IF UNSURE)"; ATMOS.PRESS 

290 INPUT "ENTER ADIABATIC FLAME TEMPERATURE IN DEGREES RANKINE"; 

FLAME . TEMP 

300 INPUT "ENTER MOLECULAR WEIGHT OF COMBUSTION PRODUCTS"; MOL.WEIGHT: 

PRINT 

310 =! eK CALCULATE VACUUM THRUST COEFFICIENT kk 

320 VACTHRUST.COEFF = CAPITAL.GAMMA * SQOR(2 / BETA) * DENOMINATOR + EXPAN- 

SION.RATIO * (1 / PRESS.RATIO) 

330 PRINT "VACUUM THRUST COEFFICIENT IS"; VACTHRUST.COEFF 

340 ! Kk CALCULATE CHARACTERISTIC EXHAUST VELOCITY KKK 

350 C.STAR = SQR(49800! * FLAME.TEMP / MOL.WEIGHT) / CAPITAL.GAMMA 

360 PRINT "CHARACTERISTIC EXHAUST VELOCITY IS"; C.STAR; "FEET PER SECOND." 

370 PRINT "PRESSURE RATIO IS"; PRESS.RATIO: PRINT 

380 INPUT "ENTER NOZZLE THROAT DIAMETER (IN INCHES)"; THROAT.DIA 

390 THROAT.AREA = 3.14159 * (THROAT.DIA / 2) * 2 

400 INPUT "ENTER THRUST LEVEL (POUNDS)"; THRUST 

410 CFPC = THRUST / THROAT.AREA 

420 CHAMBER.PRESS = (CFPC + EXPANSION.RATIO * ATMOS.PRESS) / VAC- 

THRUST. COEFF 

430 PRINT "CHAMBER PRESSURE WHEN THRUST IS"; THRUST; "POUNDS IS"; CHAM- 

BER.PRESS; "PSIA." 

1000 INPUT "DO YOU WISH TO CONTINUE (Y/N)"; ANSWERS 

1010 IF ANSWERS = "Y" OR ANSWERS = "y" THEN GOTO 10 ELSE END 
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Pyrotechnic Ignition and Propagation: A Review 


K. L. and B. J. Kosanke 
PyroLabs, Inc., 1775 Blair Road, Whitewater, CO 81527, USA 


ABSTRACT 


The ideal pyrotechnic is completely stable in 
storage and handling, yet performs its mission 
completely with absolute reliability upon de- 
mand. Many accidents in pyrotechnics are the 
result of unintentional ignitions during handling 
and storage. There can also be serious safety 
ramifications of ignition and propagation fail- 
ures of pyrotechnic devices. This review article 
presents a fairly rigorous, but mostly non- 
mathematical discussion of the ignition and 
propagation processes. 


Keywords: ignition, propagation, heat of 
reaction, activation energy, spontaneous 
ignition, thermal run-away, cook-off 


Introduction 


An understanding of the mechanism of pyro- 
technic ignition and propagation will improve 
one’s ability to identify and solve problems with 
ignition failures (duds) and unintended ignitions 
(accidents). In addition, many of these same 
principles play an important role in understand- 
ing the control of pyrotechnic burn rates. This 
article will examine these important topics 
thoroughly, however, not at a mathematically 
rigorous level. For more detailed and rigorous 
discussions, readers are referred to the writings 
of Merzhanov and Abramov.""7! 


Pyrotechnic Reaction 
Energy Considerations 


Pyrotechnic compositions are mixtures of 
fuel(s) and oxidizer(s) and possibly other materi- 
als. They are used to produce energy on de- 
mand in the form of heat, light, sound, etc. Pyro- 
technic compositions are said to be in a “meta- 
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stable” state. That is to say under typical condi- 
tions they are stable and do not react to release 
their internal chemical energy unless externally 
stimulated in some way. Probably the most 
common stimulus is the addition of heat, such 
as provided by a burning match or fuse. Igni- 
tion is the process of stimulating a pyrotechnic 
composition to release its internal energy and 
can be defined as “the initiation of self- 
sustained burning or explosion of a pyrotechnic 
material”.!! 


Figure | illustrates the process of ignition by 
graphing the internal energy of a tiny portion of 
pyrotechnic composition during the progress of 
its chemical reactions and is typical of non- 
spontaneous exothermic chemical reactions.!”! 
At the left of the graph, where the process be- 
gins, the pyrotechnic composition has a certain 
amount of internal energy. To accomplish igni- 
tion, external energy is supplied, such as from a 
burning match. This addition of energy increases 
the internal energy of the composition and is 
seen as a rise in the curve of Figure 1. This is 
indicated as the “Energy In” part of the reac- 
tion. As the process continues, eventually the 
pyrotechnic composition ignites to release its 
stored chemical energy to the surroundings. 


' Activation Energy 


Heat of Reaction 


Internal Energy 


Energy Energy 
Int 


Reaction Progress 


Figure 1. A graph illustrating the flow of 
energy into and out of a tiny portion of 
pyrotechnic composition. 
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Table 1. Heats of Reaction of Binary Pyrotechnic Compositions. 


(%) | Oxidizer 
.Magnesium (37)... 
.Magnesium (40) |. 


Magnesium 
Aluminum 


Aluminum 


This loss of internal energy is seen as a drop in 
the curve of Figure 1 and is indicated as the 
“Energy Out” part of the reaction. The energy 
that was required to stimulate this release is 
commonly referred to as the “Activation Energy” 
(E,). The net amount of energy produced by the 
pyrotechnic reaction is referred to as the “Heat 
of Reaction” (AH,).!"! 


Even the smallest particles of fuel and oxi- 
dizer in the pyrotechnic composition, are clusters 
of many billions of atoms bound together to 
form the particle. It is possible to think of the 
two-step process, energy in and energy out, as 
first when old chemical bonds are being broken 
in the fuel and oxidizer, and second, when new 
chemical bonds are being formed to make the 
reaction products. This also helps to make it 
clear why the activation energy"! requirement 
acts as a barrier that must be surmounted to ini- 
tiate the chemical reaction. Until the necessary 
energy is supplied to break the original chemi- 
cal bonds, thus freeing individual fuel and ox1- 
dizer atoms, they are not available to react with 
each other to form new chemical bonds. 


In a pyrotechnic chemical reaction, a net 
amount of energy will be produced, providing 
the new chemical bonds being formed in the 
reaction products are stronger than the old 
bonds that must first be broken in the fuel and 
oxidizer. Table | is a listing of the heats of re- 
action for some two-component pyrotechnic 
reactions. The reason that varying amounts of 
energy are produced is that in each case differ- 
ent numbers and strengths of chemical bonds are 
broken and formed. 


A collection of atoms, such as those bound 
together in a particle of fuel or oxidizer, are not 
held in absolutely rigid positions. The individ- 
ual atoms jostle about (vibrate), back and forth, 
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Iron(Il) oxide 


[5a] 


(%) AH, (kcal/g) 


Potassium chlorate (63) A298 
Potassium perchlorate (60) Jo 2.24 
Barium nitrate 

Potassium perchlorate 


and up and down. Because of these internal mo- 
tions, the individual atoms possess energy, of- 
ten referred to as thermal energy. In the process 
of jostling with one another, the atoms transfer 
some of their thermal energy from one to an- 
other. The net result of this jostling and energy 
sharing is that some atoms have much energy 
while others have little, and an atom that has 
much energy now may have little energy later. 


Figure 2 is a graph illustrating the distribu- 
tion of thermal energies! of individual atoms 
in fuel and oxidizer particles in a pyrotechnic 
composition at some temperature 7|. The curve 
goes through the origin of the graph, meaning 
that zero atoms have zero energy. Thereafter an 
increasing number of atoms have increasing 
energy, until a peak is reached, followed by 
continuously decreasing numbers of atoms pos- 
sessing higher and higher energies. Also shown 
in Figure 2 is the activation energy E,, which is 


Number of Atoms 


Thermal Energy 


Figure 2. A graph illustrating the distribution 
of thermal energies of atoms in a pyrotechnic 
composition. 
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required to initiate the pyrotechnic reaction. 
Note that the composition contains some atoms 
with energies exceeding the activation energy 
barrier. (In Figure 2 the number of atoms with 
energies greater than E, has been exaggerated 
for clarity.) 


Since some atoms in the pyrotechnic com- 
position have sufficient energy to react, the 
question should be, ““Why isn’t the pyrotechnic 
composition reacting?” The answer is, “It is 
reacting, but very, very slowly.” To see why 
this is the case, consider the following: At room 
temperature, only | atom in roughly every mil- 
lion billion (10'°) has the needed activation en- 
ergy.'*! Additionally, it is only the fuel atoms 
that are in direct contact with oxidizer atoms 
that can react. When considering the fraction of 
atoms in a tiny particle that are on its surface, 
and the fraction of surface atoms that are likely 
to be in direct contact with the right atoms on 
the surface of other particles, only 1 atom in 
roughly every thousand billion billion (107') is 
capable of reacting at any given time. 


Thermal Run-Away and 
Spontaneous Ignition 


Ifthe temperature of the pyrotechnic compo- 
sition is raised, from 7; to a higher temperature 
To, as illustrated in Figure 3, on average the at- 
oms jostle around with more energy. More sig- 
nificantly, however, the number of atoms with 
energies exceeding the activation energy barrier 
increases greatly.!"! As a consequence, there are 
now many more atoms capable of reacting, and 
there is a corresponding increase in the rate at 
which the reactions occur. Recall, however, that 
these chemical reaction produce thermal en- 
ergy; thus an increase in the reaction rate causes 
an increase in the rate of production of heat; 
which would seem to produce a further increase 
in temperature; which causes still more atoms to 
have energies exceeding the activation energy 
barrier; which causes a still greater increase in 
reaction rate and the rate of heat production; 
which causes a further increase in temperature; 
etc. This accelerating cyclic process is outlined 
in Figure 4 and leads to what can be called 
“thermal run-away” and ignition. 
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T,<T, 


Number of Atoms 


Thermal Energy 


Figure 3. An illustration of the effect of 
increasing temperature on the distribution of 
the thermal energy of atoms. 


Taken literally, the process outlined in Fig- 
ure 4, suggests that the slightest temperature 
rise of a pyrotechnic composition will eventu- 
ally lead to thermal run-away and ignition. Ob- 
viously, this is not correct. The reason is that, 
thus far, only the rate of thermal energy pro- 
duction has been considered, which is only half 
of the total picture. The rate of thermal energy 


Increase Temperature 
L-—_—_——“— 
Number of atoms with energy 
greater than E, increases 


. 
Reaction rate increases 
L 
Rate of energy production increases 
Ay 
Temperature of material rises 


L 


Repeat of cycle over and over again > 


L 


Thermal run-away = Ignition 


Figure 4. Outline of the accelerating cyclic 
process leading to thermal run-away and 
ignition. 
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loss from the pyrotechnic composition to the 
surroundings must also be considered. This 
more complete energy picture is presented in 
Figure 5, with both heat-gain and heat-loss rates 
plotted as a function of temperature. 


fx Gain 


Loss 


Rates of Heat 
Gain and Loss 


Temperature 


Figure 5. Graph illustrating the rates of heat 
gain and loss for a pyrotechnic composition as 
a function of temperature. 


The rate of the pyrotechnic chemical reac- 
tion k follows an exponential relationship,'! 
sometimes referred to as the Arrhenius equation: 


pode” (1) 


where A and R are constants, and 7 is absolute 
temperature. The rate of heat gain Q, is just the 
reaction rate multiplied by the heat of reaction: 


QO, =k * AH, (2) 


Thus, in Figure 5, the rate of heat gain curve 
passes through the origin and rises ever more 
steeply with increasing temperature. 


For a mass of pyrotechnic composition, heat 
loss from the surface will primarily be from 
convection through contact with the air. How- 
ever, any heat generated internally, will first 
need to be conducted to the surface. Accord- 
ingly, for spontaneous heat generation, tem- 
peratures at the center of the mass would nor- 
mally be highest.'! The rate of heat loss Q, from 
the center of the pyrotechnic composition de- 
pends on the thermal conductivity of the com- 
position and any packaging, the convective heat 
loss coefficient, the geometry of the sample (or 
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item), and the difference in temperature be- 

tween the center of the composition T and am- 

bient temperature 7. This may be expressed 
-[6] 

as: 


Q;= C(T— T,) (3) 


where C is a constant derived from the geome- 
try and thermal properties of the pyrotechnic 
sample (or item). Accordingly, in Figure 5, the 
rate of heat loss curve is a straight line crossing 
the temperature axis at ambient temperature and 
with a slope equal to C. 


To illustrate why, under typical storage condi- 
tions, pyrotechnic compositions are meta-stable 
and do not spontaneously ignite, as suggested 
by the process outlined in Figure 4, consider Fig- 
ures 6 and 7. Figure 6 is an enlarged view of the 
low temperature region from Figure 5. If a py- 
rotechnic composition is formulated from mate- 
rials at ambient temperature, the composition 
will be at the same temperature, at least ini- 
tially. In Figure 6, note that at ambient tempera- 
ture, the rate of heat loss is zero, while the rate 
of heat gain is greater than zero. Accordingly, 
the temperature of the sample will begin to in- 
crease. The temperature of the sample will con- 
tinue to rise until the rates of gain and loss are 
equal. This occurs at the crossing point of the 
“gain” and “loss” curves, where the temperature 
of the sample 7, has risen to slightly above ambi- 
ent temperature. (Note that in Figure 6, the 
temperature difference between 7, and T, has 
been exaggerated for clarity.)!"! 


Now imagine that for some reason the pyro- 
technic composition were momentarily raised 
from temperature 7, to 7,, somewhat further 


Loss 


Rates of Heat 
Gain and Loss 


Ta Ts 
Temperature 


Figure 6. Enlarged view of the low tempera- 
ture region of Figure 5. 
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above ambient, see Figure 7. In this case, both 
the rate of heat gain of the pyrotechnic compo- 
sition and the rate of heat loss from the compo- 
sition, increase. However, the rate of loss is 
greater than the rate of gain. Accordingly, the 
net effect will be a loss of thermal energy with 
time. Thus, the temperature of the composition 
will decrease and must continue to fall to the 
temperature where the rates of loss and gain are 
equal (at temperature 7,). From the above dis- 
cussion, it can be seen that pyrotechnic compo- 
sitions are at least meta-stable, in that any small 
addition of energy will not ultimately cause 
thermal run-away and ignition. 


Figure 7 shows that if the temperature of the 
pyrotechnic composition were to be raised mo- 
mentarily to a still higher temperature 7>, that 
the results would be quite different. In this case, 
again both the rate of heat gain and the rate of 
heat loss increase. However, this time the rate 
of gain has overtaken the rate of loss. At this 
temperature, there is a net accumulation of heat, 
producing a further increase in temperature. In 
fact, this is an ever accelerating process, be- 
cause as the temperature increases, the rate of 
gain increases much faster than the rate of loss. 
In this case, the process outlined in Figure 4 
does apply and leads to thermal run-away and 
spontaneous ignition. 


Gain Rate “Se 


Loss Rate ee 


Rates of Heat 
Gain and Loss 


Ts 


Temperature 


Figure 7. Illustration of the effect of raising the 
temperature of a pyrotechnic composition. 
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Thermal Run-Away Temperature!*! 


Obviously the temperature at which thermal 
run-away can occur for each pyrotechnic com- 
position is of great importance from a safety 
standpoint. Whenever a pyrotechnic composi- 
tion is raised above this temperature, it will be- 
gin to undergo thermal run-away and will even- 
tually ignite spontaneously. In Figure 8, the 
run-away temperature is designated as T,. and is 
the temperature corresponding to where the 
gain and loss curves cross for the second time. 
For any composition, this temperature could be 
established experimentally (with some effort 
and much time) or mathematically (providing 
the gain and loss relationships are known). 


Rates of Heat 
Gain and Loss 


T 
a Temperature 


Figure 8. Illustration of the thermal run-away 
temperature of a pyrotechnic composition, 
under one set of conditions. 


From equations | and 2 it can be seen that 
the rate of heat gain depends on the activation 
energy and the heat of reaction, which are de- 
termined by the formulation of the pyrotechnic 
composition. The heat of reaction is easy to 
calculate, providing one knows the equation for 
the chemical reaction,"! or it can be deter- 
mined experimentally.'*! Determination of acti- 
vation energy must be established through ex- 
perimentation.”'°! From equation 3, as ex- 
pressed by the constant C, it can be seen that 
the rate of heat loss depends on the thermal 
conductivity, the convection coefficient and 
geometry of the composition, and on ambient 
temperature. All of these parameters can be de- 
termined with only modest effort. However, the 
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rate of heat loss depends on many things other 
than the pyrotechnic formulation. For example, 
the degree of compaction of the composition, 
the size (mass) of the item or sample, and the 
packaging of the composition all affect the 
thermal conductivity of the composition or 
item. Also, the rate of heat loss is a function of 
ambient temperature. Thus, for each pyrotech- 
nic formulation there is not just one thermal 
run-away temperature, rather there is one for 
each of an infinite number of different condi- 
tions. This illustrates the problem in trying to 
use thermal run-away temperatures to charac- 
terize a pyrotechnic formulation (and why 
many pyrotechnists have never heard of it). 


Spontaneous Ignition Due to 
Thermal Run-Away 


Even though the use of thermal run-away 
temperature as a way of characterizing pyro- 
technic compositions is of limited value, the 
concept is important because it helps to identify 
some potentially dangerous conditions where 
there will be delayed spontaneous ignitions. For 
example, Figure 9 illustrates the effect of vary- 
ing sample size. Note that the rate of heat gain 
(per gram of composition) is unaffected by 
sample size, but the rate of heat loss is sample- 
size dependent. Small samples generally lose 
heat easily and have a rate of heat-loss curve 
that is steep, with two crossing points, the 
higher of which is the thermal run-away tem- 
perature. As the sample size increases (medium 
sample size in Figure 9), the slope of the curve 
decreases, lowering the run-away temperature 
more and more with increasing sample size. At 
some point, for a large sample, there will only 
be a single point of contact between the curves. 
This represents the largest sample, under a spe- 
cific set of conditions, that theoretically will not 
spontaneously run-away and ignite. For sam- 
ples larger than this, the rate of heat gain is al- 
ways more than the rate of loss, and the sample 
will always run-away thermally. It may take a 
very long time, but for large enough samples, 
eventually, there will always be a spontaneous 
ignition. 
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Rates of Heat 
Gain and Loss 


Very Large 


T, | Temperature 


Figure 9. An illustration of the effect of 
increasing sample size on the rate of heat 
loss and therefore on thermal run-away 
temperature. 


The rate of heat loss is sample-size depend- 
ent because sample size is one factor that af- 
fects the constant C in equation 3. Another fac- 
tor is the thermal conductivity of the pyrotech- 
nic composition and its packaging. In Figure 9, 
poorly conducting compositions and insulating 
packaging, produce effects equivalent to a large 
sample size. 


A somewhat similar situation arises for in- 
creasing ambient temperature, see Figure 10. 
When the ambient temperature is low, samples 
lose heat to the surroundings fairly easily. This 
places the heat-loss rate curve fairly high on the 
heat-gain curve, producing two crossing points, 
the higher of which is the thermal run-away 
temperature. As the ambient temperature in- 
creases, it becomes more difficult for the sam- 
ple to lose heat. The slope of the heat-loss curve 
is unchanged, but its position relative to the 
heat-gain curve is lower. This lowers the run- 
away temperature more and more with increas- 
ing ambient temperature. At some point, for 
high enough ambient temperature, there will 
only be a single point of contact between the 
curves. This represents the highest ambient 
temperature, for this type and size of sample, 
that will not spontaneously run-away thermally 
(ignite). For ambient temperatures greater than 
this, the rate of heat gain is always more than 
the rate of heat loss, and the sample will always 
run-away thermally. It may take a long time, 
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but for such hot ambient conditions there will 
eventually be a spontaneous ignition. 


Rates of Heat 
Gain and Loss 


Ta, !ay'a, la, Temperature 


Figure 10. An illustration of the effect of 
increasing ambient temperature on the rate 

of heat loss and therefore on thermal run-away 
temperature. 


Just how large a sample and just how high 
an ambient temperature are required for thermal 
run-away and spontaneous ignition depends on 
the chemical formulation and the conditions of 
its packaging and storage. For some composi- 
tions under favorable conditions, it may require 
millions of tons of material and take years to 
run-away and spontaneously ignite. However, 
under more extreme conditions, or for other 
compositions, tiny samples of composition may 
ignite very quickly. 


Time to Ignition, Cook-Off Tests 


The time for any given pyrotechnic compo- 
sition to ignite is a function of the temperature 
to which it is exposed, as illustrated in Figure 11. 
If a sample is placed in an oven, its temperature 
will begin to rise from ambient temperature T,,, 
eventually reaching the temperature of the 
oven. If the temperature of the oven T; is less 
than the run-away temperature 7, for the pyro- 
technic composition, the sample will never ig- 
nite (1.e., the time to ignition is infinite). 
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—— Ignition —— 


Temperature 


Figure 11. An illustration of the dependence of 
time to ignition on the temperature to which a 
pyrotechnic composition is heated. 


If the temperature of the oven 7) is slightly 
higher than the run-away temperature, upon 
placement in the oven, the sample will start to 
heat up and will eventually rise to the tempera- 
ture of the oven. For a while, it may appear that 
nothing is happening with the sample. However, 
inside the sample, heat is slowly accumulating, 
raising the temperature, at first perhaps imper- 
ceptibly. As the internal temperature rises, the 
rate of reaction increases, increasing the rate of 
heat gain, and further increasing the tempera- 
ture (i.e., thermal run-away has begun). As a 
result of this ever accelerating process, the inter- 
nal temperature rises ever more rapidly, until 
eventually there is an ignition, at time f, in Fig- 
ure 11. 


Ifthe temperature of the oven (73) is signifi- 
cantly higher than the run-away temperature, 
upon placement in the oven, the sample again 
will heat up, approaching the oven temperature. 
However, under these conditions, the sample’s 
temperature rise may not slow significantly as it 
reaches the oven temperature, before thermal 
run-away is at an advanced state with ignition 
occurring more quickly at time ¢3. 


Determining the time to ignition as a func- 
tion of temperature, has important ramifications 
for the storage of pyrotechnics (and explosives). 
If one were to guess wrong, the consequences 
of an accidental spontaneous ignition could be 
disastrous. Tests performed to discover the time- 
to-ignition and temperature relationship are 
sometimes called “cook-off” tests.'''! In these 
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tests, samples or items are typically placed in a 
heated bath, after having thermocouples in- 
stalled internally. Bath and sample temperatures 
are monitored as a function of time from the 
start of the test, and the time to ignition or ex- 
plosion (if either occurs) is recorded. 


Ignition and Ignition Temperature 


Ignition is one of the more difficult terms to 
define in pyrotechnics. Obviously it cannot be 
defined as when chemical reactions start. As 
was discussed earlier, some pyrotechnic reac- 
tions are occurring all the time, although at a 
very, very low rate. Even after the thermal run- 
away temperature has been reached, there may 
be no obvious sign anything is happening. For 
most observers, the appearance of a flame (high 
temperature radiant gases) or at least obvious 
incandescence of the solid phase is taken as the 
indication that ignition has occurred. As sug- 
gested in Figure 11, at the time of ignition very 
rapidly accelerating reaction rates produce a 
near instantaneous rise in temperature, typically 
from several hundred to two thousand degrees 
Celsius or higher. Thus the physical manifesta- 
tions of ignition develop very rapidly as igni- 
tion is occurring. 


Ignition temperature can be defined as “the 
minimum temperature required for the initiation 
of a self-propagating reaction”.!'”! However, 
from the above discussion, that temperature can 
vary widely depending on sample conditions 
and on how long one is willing to wait for the 
ignition to occur. These problems are mostly 
eliminated for ignition temperature measure- 
ments, because the conditions and delay time 
are usually specified in the procedure to be 
used. Unfortunately, there are many different 
procedures that are used; the Encyclopedia of 
Explosives''*! alone lists 14 different methods. 
This means there can be more than one ignition 
temperature reported for the same pyrotechnic 
composition, depending on which method was 
used. Fortunately, the various ignition tempera- 
tures of the most commonly used methods all 
tend to be in the same general range, primarily 
because the measurement conditions of the 
various methods tend to be somewhat similar. 
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Obviously, however, the most consistent results 
will be achieved if all measurements are made 
using the same method. Also, when evaluating 
reported ignition temperature data, it is useful 
to know which method has been used. Three 
hot bath methods and one differential thermal 
analysis method are described in the following 
paragraphs. 


Hot Bath Method One: Ignition temperature 
is the lowest temperature of a bath of Wood’s 
metal,'! to within 5 °C, that results in ignition 
within 5 minutes for a 0.1 g sample in a pre- 
heated small glass test tube inserted 1/3 its 
length into the bath." 


Hot Bath Method Two: Ignition temperature 
is the temperature of a bath of Wood’s metal 
that results in ignition in 5 seconds (determined 
graphically using time to ignition data) for a 
1.0 g sample in a thin-walled brass or copper 
tube (typically a No. 6 detonator shell).!"°! 


Hot Bath Method Three: Ignition tempera- 
ture is the temperature of a bath of Wood’s 
metal, heated at a rate of 5 °C per minute, at 
which ignition occurs for a 0.5 g sample in a 
tightly corked glass test tube (125 mm long by 
15 mm inside diameter with a 0.5 mm wall 
thickness).!"*! 


Differential Thermal Analysis Method: \gni- 
tion temperature is the temperature of onset of 
the ignition exotherm for a 10 to 100 mg sam- 
ple heated at a rate of 50 °C per minute.!"*! 


A collection of ignition temperatures for a 
series of two-component, stoichiometric pyro- 
technic compositions is presented in Table 2. 
Shidlovskiy reports the method as one using an 
electric furnace instead of a bath of Wood’s 
metal. Unfortunately, he provides no informa- 
tion on the method, other than an estimate of 
the accuracy of the results to be within 10 °C. 
(This data was chosen for inclusion because it 
is the most systematically complete set of data 
known to the authors.) 


It might be of interest to note that a typical 
pyrotechnic composition raised to its ignition 
temperature will have about 30 million times 
more atoms with energies exceeding E, than at 
room temperature. |*! 
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Table 2. Ignition Temperatures for Binary Pyrotechnic Compositions. 


ium chlorate 


[Se] 


Ignition Temperature (°C) 


Al dust 


The Effect of Melting and 
Tammann Temperature 


For a pyrotechnic reaction to occur, the at- 
oms (or molecules) must have the required acti- 
vation energy, and they must be in direct con- 
tact with another atom of the correct type. Even 
for well-mixed solid particles, there are rela- 
tively few points of contact between individual 
particles; see the top illustration of Figure 12. 
Thus the number of fuel and oxidizer atoms that 
are in contact with one another is normally 
quite small. However, if one of the components 
melts to flow around the surfaces of the other 
particles, there is a great increase in the number 
of atoms in contact; see the bottom illustration 
of Figure 12. 


Accordingly, melting can have a significant 
effect on the likelihood of ignition. The poten- 
tial effect on ignition is outlined in Figure 13. If 
the percentage of atoms in physical contact in- 
creases upon the melting of one component, 
then more atoms with energies exceeding the 
activation-energy barrier will be in contact with 
one another. That means the reaction rate will 
then be greater, and with it the rate of produc- 
tion of thermal energy, which means that ther- 
mal run-away and ignition can occur at a lower 
temperature. Thus it is suggested that if a com- 
position is nearing its ignition temperature, and 
one component of the composition melts, that 
could result in ignition occurring at that lower 
temperature. For example, the melting point of 
potassium nitrate and the ignition temperature of 
Black Powder are effectively the same.!'””! 


Melting can be thought of as occurring when 
the thermal vibrations of a solid are so strong 
that some of the bonds, which had been holding 
the solid together, are broken. While for most 
pure chemicals melting has a sudden onset at a 
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Before: 


Very little 
surface in 
contact 


After: 


Very much 
more surface 
in contact 


Figure 12. An illustration of the great increase 
in contact between fuel and oxidizer after one 
component melts. 


specific temperature, the vibrations of the atoms 
in the solid become increasingly strong as the 
temperature is increased toward the melting 
point. This can be thought of as the loosening 
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Melting of one component 
L 


The percentage of atoms in 
physical contact increases 


J 


More atoms with energies exceeding 
E, are now in contact and will react 


1 
Rate of reaction increases 
1 
Rate of energy production increases 
L 


Thermal run-away (ignition) occurs 
at a lower temperature 


Figure 13. An outline of the effect of melting 
on thermal run-away and ignition temperature. 


of the bonds holding the solid together. It has 
been suggested that at an absolute temperature 
halfway to the melting point, the bonds become 
so loose that there can be a significant commin- 
gling of the atoms of particles that are in direct 
contact.'**! This temperature is called the Tam- 
mann temperature after the researcher making 
this observation. 


Tammann temperatures are of interest be- 
cause samples at or above these temperatures 
demonstrate significantly increased sensitivity 
to accidental ignition. Table 3 lists the Tam- 
mann temperatures for some common pyro- 
technic oxidizers. 


Table 3. Tammann Temperatures for 
Common Pyrotechnic Oxidizers.!"*" 


Tammann 
Oxidizer Temp. (°C 


Sodium nitrate 


Potassium perchlorate 
Lead chromate 
lron(IIl) oxide 
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Propagation and the 
Propagation Inequality 


Having successfully ignited a pyrotechnic 
composition is no guarantee that the reaction 
will propagate throughout. This is because the 
application of an external stimulus, such as a 
flame, typically provides thermal energy to only 
a small portion of the composition, and the igni- 
tion stimulus is usually of relatively short dura- 
tion. After its application, the pyrotechnic com- 
bustion reaction will continue to propagate 
through the composition only so long as the 
pyrotechnic reaction provides sufficient energy 
to the unreacted composition. What is needed is 
sufficient energy to raise the unreacted compo- 
sition above its ignition temperature. This proc- 
ess is illustrated in Figure 14. The portion of the 
rod of pyrotechnic composition to the extreme 
right has already been consumed by burning. 
Just to the left of that is shown a thin disk of 
composition that has ignited and is still reacting 
(burning). Just left of that is another thin disk of 
composition, labeled “pre-reacting material”, 
which has not yet ignited. This disk of pre- 
reacting material will only ignite if it is raised 
above its ignition temperature, which means 
that a significant number of its atoms and mole- 
cules will have received at least the required 
activation energy E,. 


Unreacted 
Composition 


Consumed 


wh 


SS Reacting Material 
Pre-Reacting Material 
Figure 14. Illustration of a propagating rod of 
pyrotechnic composition. 


Most of the energy being produced by the 
reacting material (the heat of reaction AH,) is 
lost to the surroundings. However, some frac- 
tion Fy, of the energy will be fed back from the 
reacting layer to the pre-reacting layer. The ac- 
tual amount of energy fed back is just the prod- 
uct of the heat of reaction and this fraction (i.e., 
AH,*F »). Propagation will occur providing more 
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energy is fed back than is required for ignition 
of the pre-reacting disk of composition. This 
statement, then, gives rise to something that 
could be called the “propagation inequality” 


AH, * Fy > Ei (4) 


Obviously, the probability of successful py- 
rotechnic propagation is increased by anything 
that results in more energy being produced by 
the burning composition, a greater percentage 
of that energy being fed back, or a reduction of 
the activation energy requirement. The heat of 
reaction and the activation energy are deter- 
mined by the chemical nature of the composi- 
tion; however, a thorough discussion is beyond 
the scope of this article.!*! Energy is fed back to 
the pre-reacting layer through any combination 
of conduction, convection, or radiation. These 
are more fully described in Figure 15, where 
some of the factors influencing the efficiency of 
these mechanisms are also presented. 


Figure 16 is similar to Figure 14 but pro- 
vides a more compete description of the propa- 
gation process, including information on the 
relative temperatures expected. Zone (a) is un- 
reacted pyrotechnic composition, which has thus 
far been unaffected and remains at ambient 
temperature. Zone (b) is described as the warm- 
up zone, where the temperature has started to rise 
above ambient, as a result of thermal conduc- 
tion and possibly convection. It is in this zone 
where reaction rates are first beginning to in- 
crease. These reactions are sometimes referred to 
as “pre-ignition reactions”**"! and contribute 
relatively little thermal energy. In Zone (ce), the 
temperature has risen significantly, at least one 
component of the composition has melted, and 
some gaseous materials may be bubbling to the 
surface. Because of the rise in temperature and 
the greatly increased contact between fuel and 
oxidizer, the reaction rate in Zone (c) and the 
production of heat is greatly increased. In 
Zone (d), much of the reaction is occurring in 
the gas phase; however, some droplets of react- 
ing composition ejected from the surface may 
be present. Again the reaction rate and the 
thermal energy being produced has increased 
substantially from the previous zone, and the 
temperature has peaked. In the final region, 
Zone (e), the energy producing reactions have 
ceased, and as a result of heat loss to the sur- 
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Conduction: 


— Thermal energy (molecular 
vibration) is conducted along sol- 
ids, from hotter to cooler. 


— Factors maximizing conductive 
feedback: 
¢ Compacted composition 
¢ Metal fuels 
¢ Metal casing or core wire 


Convection: 


— Hot gases penetrate the solid 
composition along spaces 
between grains (fire paths). 


— Factors maximizing convective 
feedback: 


¢ Uncompacted composition 
¢ Granulated composition 


Radiation: 


— Thermal (infrared) radiation, 
emitted from flame and glowing 
particles, is absorbed by incom- 
pletely reacted composition. 


— Factors favoring radiative 
feedback: 


¢ Solid or liquid particles in flame 


¢ Dark or black composition 


Figure 15. Outline describing the mechanisms 
of pyrotechnic energy feed back and the factors 
that affect them. 


roundings, the temperature begins to fall sig- 
nificantly. 
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Flame 
Temperature 


Ambient 


——_—_ 


Temperature 


Distance Along Composition 


Figure 16. A more complete description of the 
processes of pyrotechnic propagation: a, 
unreacted composition; b, warm-up zone; ¢, 
condensed phase reactions; d, gas phase 
reactions; and e, reactions complete. 
(Reference based on 5d.) 


Conclusion 


The subjects of ignition and propagation 
could have been dealt with much more expedi- 
tiously. Specifically, it could simply have been 
stated that ignition requires raising at least a 
portion of a pyrotechnic composition to its igni- 
tion temperature, and that propagation requires 
the feed back of sufficient energy for continuing 
self-ignition of the composition. However, while 
this would have saved much time, it would have 
provided little understanding of the important 
principles involved. Accidents caused by unin- 
tentional ignitions continue to plague the pyro- 
technics industry. In addition, ignition failures 
with fireworks, whether leaving an unfired-aerial 
shell in a mortar that needs to be cleared or re- 
sulting in a dud shell left after a fireworks dis- 
play, have serious potential for accidents as well. 
Ignition failures with signaling smokes or flares 
may cause someone to not be rescued. It is 
through a more thorough understanding of the 
basis for ignition and propagation, that pyrotech- 
nists will be better able to solve potential prob- 
lems before accidents occur. 
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Notes 


[a] More recent usage of this term is “Enthalpy 
of Reaction”. 


[b] The use of the term “activation energy” in 
the context of solid-phase pyrotechnic com- 
positions has a slightly different meaning 
than in aqueous or gas-phase chemistry. In 
aqueous and gas-phase chemistry, molecules 
can be thought of as reacting individually, 
or at least in small groups. In that case, ac- 
tivation energy can be thought of as the 
amount of energy needed for a collision be- 
tween two individual molecules to cause 
them to react. However, a typical pyrotech- 
nic composition is composed of solid parti- 
cles of fuel and oxidizer, with each particle 
comprised of many billions of atoms or 
molecules. In this case, activation energy 
takes on much more of a macroscopic and 
less precise meaning. For the purpose of this 
article, activation energy of solid-phase py- 
rotechnic compositions is generally taken to 
mean that amount of thermal energy needed 
to induce a sustained exothermic reaction 
within a tiny portion of the composition. 


[c] Pyrotechnic compositions, especially those 
of dark color, exposed to bright sunlight 
may be an exception to the generalization 
about the highest temperature expected to 
be in the middle of the composition. 


[d] Note that it is assumed the composition is 
physically at rest (i.e., it is not being mixed, 
which adds energy to the system). 


[e] In some texts thermal run-away temperature 
is called the critical temperature!" or the re- 
action temperature.'*"! 


[f] Wood’s metal is a eutectic alloy of bismuth, 
lead, tin and cadmium. It melts at 70 °C. 


[g] Some oxidizers are known for producing 
compositions with low activation energies 
(e.g., chlorates), while other oxidizers tend 
to produce compositions with high activa- 
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pa) 


tion energies (e.g., oxides and sulfates). Fu- 
els with low melting points or low decom- 
position temperatures (e.g., sulfur, lactose 
and acaroid resin) tend to form compositions 
with low activation energies. Metal fuels 
tend to produce compositions with high 
heats of reaction. 


Pre-ignition reactions are typically reactions 
taking place in the solid state. While such 
reactions can be a source of energy, they 
generally only contribute in a minor way to 
promoting ignition. This is because solid- 
state reaction rates are constrained by the 
difficulty of fuel and oxidizer commingling 
while both remain solid. 
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ABSTRACT 


The presence of sulfur/chlorate mixtures in 
fireworks compositions has been proposed as 
the cause of numerous accidents in production 
and storage. A series of accidents after the in- 
troduction of the 1875 Explosives Act led to 
additional legislation to prohibit such mixtures 
in the UK. In this paper, as part of an ongoing 
programme to quantify the hazards associated 
with the presence of the mixture in fireworks 
compositions, we have reviewed previously 
published work and give an outline of the ap- 
proach we will take in our current research. 


Keywords: sulfur, chlorate, pyrotechnics, 
fireworks, reactivity, accidents, safety 


Introduction 


Sulfur/chlorate mixtures have, for many 
years, been known as potential causes of the 
accidental initiation of pyrotechnics. In the pe- 
riod 1875 to 1894 there are numerous incidents 
recorded in the United Kingdom (UK) HM In- 
spectors of Explosives reports that are attrib- 
uted to the presence of sulfur/chlorate mixtures. 
The majority of these incidents, which resulted in 
11 deaths along with many injuries, were be- 
lieved to be caused by the spontaneous ignition 
of the mixture while the remainder were thought 
to be due to the sensitiveness''! of the mixtures. 


An accident"! in 1887 was attributed to a 
green star composition where both barium and 
potassium chlorates were present with sulfur 
(KCIO;3 47%, Ba(ClO3)2 9%, Ba(NO3)2 24%, 
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S 14%, shellac 4%, charcoal 2%), whereas an 
accident’! in 1890 was attributed to a green 
lance manufactured using two compositions. 
These comprised a green fire which contained 
both potassium and barium chlorates and a 
white fire containing sulfur. It was believed that 
filling was not carried out carefully and that the 
two compositions were able to mix and form a 
sulfur/chlorate mixture. Similar ignitions were 
reported"! involving a green fire composition 
containing 62% barium chlorate with 8% sulfur. 


The acidic nature of the materials used at 
that time was proposed as a possible cause of 
the accidents.) The occurrence of such events 
led to an Order in Council in 1894 prohibiting 
the manufacture, import, keeping, conveyance 
or sale of fireworks containing mixtures of sulfur 
with any chlorate except with the consent of a 
Government Inspector. Such restrictions on 
fireworks have been and continue to be in place 
for fireworks manufactured in or imported into 
the UK. 


In a more recent review of accidents, in this 
case to minors from chemical experimentation 
in the 1950s, Burns'®! reported that 36% of the 
incidents were attributed to handling sul- 
fur/chlorate mixtures in conjunction with car- 
bon. 


In recent years a number of fireworks have 
been examined by the Health and Safety Labo- 
ratory for the presence of sulfur/chlorate mix- 
tures, both qualitatively’! and quantitatively.'*! 
A proportion of the fireworks tested was found 
to contain sulfur/chlorate mixtures to different 
extents within a component composition, and in 
many cases this necessitated both removal of 
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the item from sale and destruction of the fire- 
works. While there have been no recent re- 
ported manufacturing or storage incidents in the 
UK, there are reports of accidents with pyro- 
technic materials in many countries (e.g., 
Peru,”! China!'®!! and India!’?! that could be 
attributable to the presence of sulfur/chlorate 
mixtures.) 


Role of Sulfur and Chlorate in Fireworks 


Sulfur is one of the commonly used fuels in 
fireworks manufacture and is said to act as a 
“tinder”, and its low temperature exothermic 
reactions with a variety of oxidisers are used to 
initiate other higher temperature reactions. '"*! 


Chlorates are present in compositions to ei- 
ther: 


a) act as an oxidiser and/or 
b) act as a chlorine source. 


The presence of chlorate is said to generate com- 
positions with relatively low ignition tempera- 
tures, in the region of 200 °C for many compo- 
sitions.''*! Conkling reports a strong exotherm 
for sulfur/potasstum chlorate mixtures corre- 
sponding to ignition at temperatures well below 
200 °C.'"*! It is further suggested by Conkling 
that the low Tammann temperature (defined as 
0.5 T where T is the melting temperature in 
Kelvin) of 42 °C, the exothermic decomposition 
of potassium chlorate, and the melting tempera- 
ture of sulfur of 119 °C are responsible for the 
low ignition temperature. At the Tammann 
temperature, a solid has 70% of the vibrational 
freedom present at the melting temperature and 
considerable diffusion into the lattice is likely. 
Similar ignitions have been observed with thio- 
urea/potassium chlorate mixtures in process- 
ing!'®! where reactivity was attributed to migra- 
tion of loosely bound sulfur from thiourea into 
the potassium chlorate lattice producing a mate- 
rial that was both highly sensitive to mechanical 
stimulus and had temperatures of ignition of 
about 150 °C. 
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Previous Work 


Several previous studies have investigated 
the reaction of sulfur and metallic chlorates. All 
workers have found the mixtures to be readily 
ignited at low temperature or with low me- 
chanical stimulus. As early as 1881, acidity of 
the sulfur was cited by Dupré!’ as a major con- 
tributing factor in the auto-ignition of sul- 
fur/chlorate mixtures. Additionally, where bar- 
ium salts were present under moist conditions 
with potassium chlorate, a double decomposi- 
tion reaction could occur to produce barium 
chlorate. Barium chlorate was also shown to be 
more susceptible to auto-ignition than the po- 
tassium salt when mixed with sulfur. Later 
work by Dupré!'*! concluded that under condi- 
tions of heat retention it was a matter of when 
not if mixtures containing sulfur and chlorate 
would spontaneously ignite. 


Amiel''”! attributed the spontaneous ignition 


of chlorates and sulfur to the production of 
chlorine dioxide (C102). Chlorates, other than 
alkali-metal chlorates, were shown to produce 
ClO, at temperatures around 60 °C in the pres- 
ence of water. Additionally, the chlorates of 
magnesium, nickel, cobalt, cadmium and zinc 
were reported to react very suddenly. 


Taradoire’”! focused his studies on potas- 
sium, barium and lead chlorates. Under wet 
conditions the production of sulfuric acid fol- 
lowed by the formation of chloric acid (HCIO3) 
was proposed as the cause of spontaneous igni- 
tion of mixtures. The mixtures formed by sulfur 
with the chlorates of barium and lead were said 
to ignite at “ordinary” temperatures. 


Ina further paper, Amiel"! reported the ef- 
fect of sample size, indicating that for large 
(unspecified) quantities, self-reaction was im- 
mediate rather than taking minutes or hours as 
for smaller (unspecified) quantities. 


Tanner'”! has demonstrated that sulfur diox- 
ide can initiate sulfur/chlorate mixtures. He pro- 
posed that the exothermic reaction: 


SO,+ 2 KC1O;+4S > 
2 SO, + S2Cl + K,SO4 


was responsible for the ignition of the material. 
This chain reaction would be capable of self 
sustaining once the initial sulfur dioxide had 
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been generated from polythionic acids (H2S,O4) 
by the equation: 


H2S,0¢6 > H,SO, a SO, “1 (n—2) S 


Polythionic acids form on the surface of sul- 
fur crystals, particularly at elevated tempera- 
tures due to surface oxidation. Decomposition 
of polythionic acids occurs when evaporation of 
surface moisture or a temperature rise make the 
acid unstable. 


An alternative mechanism was proposed by 
Rudloff"?! from thermal analysis work carried 
out on the reaction of sulfur vapour and potas- 
sium chlorate. The overall reaction for the com- 
plex mechanism was reported as: 


3 KCIO;+3S—> 
K,SO,4 + KCI + 2 SO2 + CLO 


Storey"! reported friction sensitiveness of 


sulfur/chlorate mixtures and the results of a 
number of thermal techniques that were used to 
assess the likelihood of the ignition of sul- 
fur/chlorate mixtures. Stoichiometric mixtures 
of sulfur and potassium chlorate were shown to 
have a Figure of Friction 0.12, compared to 
0.84 for lead azide. Since low Figure of Friction 
values indicate a more sensitive composition, 
this illustrates that under certain circumstances 
sulfur/chlorate mixtures can be more sensitive 
to friction than primary explosives commonly 
used as initiator compositions. Mixtures contain- 
ing less than 2.5% sulfur with potassium chlo- 
rate were shown to be insensitive to friction. 


Hot stage microscopy of mixtures contain- 
ing 1—5% sulfur with potassium chlorate indi- 
cated that they started to decompose at tempera- 
tures in the range 165-178 °C, dependent on 
sample size and percentage of sulfur. Differen- 
tial thermal analysis on 20 mg samples reduced 
the initiation temperature to the range 130-155 
°C for similar mixtures, again depending on the 
percentage of sulfur and, in this case, the source 
of the sulfur. Times to autoignition for 
stoichiometric mixtures were also considered. 
Sulfur/potassium chlorate mixtures (1 g sam- 
ples) were sealed in glass test tubes and the 
time to ignition was recorded. Sulfur from two 
sources produced similar results taking 2.5—3 
hours for ignition at 90 °C, while the sulfur 
from a third source ignited in 13.5 minutes at 
the same temperature. The sulfur samples were 
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from three different UK manufacturers and were 
representative of the materials used in fireworks 
manufacture. Unfortunately no measurement of 
acidity was reported. 


Calorimetry (Seteram C80) on 0.5 g samples 
of stoichiometric sulfur/potassium chlorate mix- 
tures indicated that exothermic reactions were 
initiated within a few hours at temperatures be- 
low 60 °C under isothermal conditions. Using a 
heating rate of 0.1 °C min !, exotherms were 
observed for similar mixtures at temperatures as 
low as 62 °C. Mixtures with low sulfur percent- 
age showed higher onset temperatures, but gen- 
erally these were below the sulfur melting tem- 
perature. The addition of sulfur dioxide to a 
mixture was shown to initiate the reaction al- 
most instantaneously. 


Research Plan 


Thermal Studies 


The early papers by Amiel and Tara- 
doire!!’”°*"! do not report the sizes of the sam- 
ples taken, the ratio of sulfur to chlorate in the 
experiments, nor the conditions under which the 
experiments were carried out. However, they do 
indicate an early appreciation of the thermal 
instability of this mixture. Tanner’s work!” util- 
ised an almost stoichiometric mixture (33:67) and 
the acidity of the sulfur was measured. 


A previous study,”“! which defines the con- 
ditions employed for the experimental work, 
used samples contained in glass or metal con- 
tainers, many of which were sealed systems. 
The use of a sealed system will prevent the dis- 
persion of any sulfur dioxide produced which 
could in turn accelerate the autocatalytic reac- 
tion yielding unrealistically low reaction tem- 
peratures or times to ignition. Tanner’s work,|”7! 
on the other hand, utilised cardboard tubes 
made from match books after the striker com- 
position had been tor off. This provided an 
open containment system and the samples tested 
would therefore not suffer from any effects of 
confinement. However, his use of a heat lamp 
and thermocouple positioned in the centre of 
the sample may have led to unrepresentatively 
low measured temperatures of ignition since the 
reaction surface of the material would have 
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been at a significantly higher temperature than 
that indicated by the thermocouple. Addition- 
ally, it may be possible that some phosphorus 
contamination remained from the matches. 


Our studies will investigate the initiation of 
sulfur/chlorate mixtures in containers that are 
similar to those used in the fireworks industry 
(i.e., cardboard tubes sealed with a clay plug at 
the bottom and either open at the top or sealed 
with a further clay or cardboard disk.) Tempera- 
ture measurement of the samples will be made 
using a thermocouple inserted into the composi- 
tion. The cardboard tubes will be heated in a 
metal block to achieve uniform temperature 
throughout the sample. To fully understand the 
reactivity of sulfur/chlorate mixtures in fire- 
works a number of variables will be systemati- 
cally examined. These include: 


1. the acidity of the sulfur, 

. the stoichiometry of the mixture, 

. the size of sample, 

. the degree of sample compaction, 

. the moisture content of the mixture, 


Nn BW WN 


. the particle size of the components of the 
mixture, and 


7. the inclusion of other fireworks materials 
such as barium or strontium salts and 
metal particles. 


Sensitiveness Studies 


While the majority of accidental ignitions of 
compositions containing sulfur/chlorate mix- 
tures have been attributed to thermal instability, 
Dupré”! reported that he could initiate a sul- 
fur/chlorate containing “green Roman” com- 
position by friction using the wooden plunger 
used in filling, particularly when the material 
had been heated. He concluded that the compo- 
sition “possesses a degree of sensitiveness 
which renders it in the highest degree a hazard- 
ous composition to deal with especially in hot 
weather”. Therefore, in parallel with our pro- 
posed thermal studies, the sensitiveness of dif- 
ferent sulfur/chlorate mixtures to mechanical 
stimuli will also be quantified using the same 
matrix of variables identified for the thermal 
stability work. 


Friction has previously been shown to be a 
major cause of accidental ignitions for explo- 
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sives in general””*! and pyrotechnics!” in par- 


ticular. However, for completeness both friction 
and impact sensitiveness will be examined. This 
will be carried out using the BAM friction and 
BAM Fallhammer apparatus to gain data using 
the UN recommended test methods” for fric- 
tion and impact, respectively. 


Future Work 


This paper has outlined the problems associ- 
ated with the presence of sulfur/chlorate mix- 
tures in fireworks and forms the introduction to 
a series of papers where we intend to report the 
findings of our research programme. Part 2!*! 
will deal with the early work carried out on dry 
stoichiometric mixtures. Future papers will 
cover other variables considered in the study. 
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ABSTRACT 


The use of chlorine-free high performance 
oxidizers in solid rocket propellants may reduce 
launch costs, while producing less polluting 
exhaust products. Hydrazinium nitroformate 
(HNF) is a candidate for such a next generation 
oxidizer. In this article the properties of HNF 
are addressed. Its combustion and decomposi- 
tion characteristics, and its application in 
rocket propellants are discussed. 


Keywords: oxidizer, hydrazinium 
nitroformate, HNF, monopropellant, 
composite propellant, decomposition 


Introduction 


Solid rocket propellants can be divided into 
two main categories. Double-base propellants, 
which are homogeneous mixtures, with main 
ingredients nitroglycerin and nitrocellulose. 
Another group consists of the so-called compos- 
ite propellants, which are heterogeneous mix- 
tures of oxidizer and fuel. The oxidizer is an 
inorganic salt, which is held together by a 
polymer matrix, which also acts as a fuel. Dou- 
ble-base propellants are mostly used in small 
military missiles. The mechanical properties 
and performance of composite propellants are 
better than those of double-base propellants. 
Therefore, most larger systems (e.g., the boost- 
ers of the Space Shuttle) are based on compos- 
ite propellants. 


The majority of the composite propellants 
use hydroxyl-terminated polybutadiene (HTPB, 
~CH,5) as binder and ammonium perchlorate 
(AP, NH,ClO,) as the oxidizer.''! Often fine 
aluminum is added as a fuel to increase the 


Page 36 


combustion temperature and hence perform- 
ance. These types of composites are cheap, reli- 
able, and can be tailored to a large range of re- 
gression rates by using burn-rate modifiers and 
varying AP particle-size distribution. However, 
there is a lot of research on new oxidizers to 
replace AP. Reasons for this are: 


¢ Byreplacing AP composites with more en- 
ergetic composite propellants, launch costs 
can be decreased dramatically. Reduction 
of launch costs is important in the current 
period of increased competition between 
launcher organizations. 


e AP-based composites are characterized by 
a large amount of hydrochloric acid (HCl) 
in their exhaust products. During each 
launch of the Space Shuttle, 162 metric 
tons of HCl are ejected into the earth’s at- 
mosphere! HCl is known for it ozone de- 
pletion properties. 


¢ The presence of HCI causes a large smoke 
trail behind missiles. This so-called secon- 
dary smoke consists of condensed water. 
(Primary smoke is caused by solid parti- 
cles in the exhaust products, for example, 
aluminum oxide and soot.) In the presence 
of HCl, water condenses more easily, as 
HCl decreases the effective saturation 
pressure of water (like the white trail be- 
hind airplanes). The secondary smoke in- 
creases the signature of a missile and is 
therefore unwanted. 
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New Oxidizers 


The introduction explains the search for new 
oxidizers that may replace AP. Several alterna- 
tives have been identified, and some of them 
have already been applied in actual systems. 


Ammonium nitrate (AN, NH,NO;) is an oxi- 
dizer that does not contain chlorine. AN-based 
composite propellants are characterized by their 
low flame temperature. This makes them ideal 
for gas generator applications. However, their 
performance is lower than that of AP-based 
composite propellants. 


The cyclic nitramines, cyclotrimethylenetri- 
nitramine (RDX, C3HsN,6Os) and cyclotetra- 
methylene tetranitramine (HMX, C,HgN3Os), 
are not real oxidizers, as they are stoichiometri- 
cally balanced. Both nitramines have a high 
flame temperature and generate combustion 
products with a low molecular weight. By re- 
placing AP in AP-based composites with these 
nitramines, the performance may be increased. 
RDX- and HMX-based propellants do not show 
a strong particle-size dependent burning rate.!”! 
This has been attributed to the melt layer formed 
during decomposition of these materials. Re- 
gression rates of nitramine-based composites 
can therefore not be tailored by particle-size 
distribution control. Also combustion catalysts 
have only been proven partially successful. AP 
is therefore often only partially replaced by 
HMX or RDX. The performance gain of nitra- 
mine composites is still not the maximum that 
can be obtained as the binder decomposition 
products cannot be oxidized, because the nitra- 
mines are stoichiometrically balanced. 


A new energetic filler is CL20 or hexanitro- 
hexaazaisowurtzitane (HNIW, CeH6Ni20)2). 


Table 1. Properties of Several Oxidizers and Energetic Fillers. 


Unlike RDX and HMX, CL20 is not yet used in 
operational systems. Its main advantage over 
the cyclic nitramines is its high density, 2.04 
g/em*. CL20 research is carried out in France 
and the United States. 


Ammonium dinitramine [ADN, 
NH,N(NO;),] and hydrazinium nitroformate 
[HNF, N2H,C(NO;)3] are two “new” high- 
performance oxidizers that do not contain chlo- 
rine and have a positive oxygen balance. In the 
United States, research focuses on the devel- 
opment of ADN-based propellants. There are 
also ADN programs in France and Sweden. 
There are indications that Russian rocket en- 
gines based on ADN have already been opera- 
tional for many years."! However, ADN is 
hampered by its low melting point (93 °C) and 
hygroscopicity."! 


HNF was part of a large research project for 
new high-performance propellants in the 1960s 
in the United States. HNF appeared to be in- 
compatible with the double bonds in unsatu- 
rated binders like HTPB."! The recent commer- 
cial production of saturated energetic binders, 
like glycidyl azide polymer (GAP), led to a re- 
newed interest in HNF. The theoretical per- 
formance of HNF-based composites is un- 
matched by any other energetic material. Be- 
cause of HNF’s high performance, the European 
Space Agency (ESA) decided to fund HNF re- 
search. 


Several properties of the oxidizers discussed 
are summarized in Table 1.'°”! The performance 
of some typical propellants with some of these 
oxidizers is given in Figure 1. 


[6,7] 
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Figure 1. Performance of several propellants 
(Chamber pressure 10 MPa, expansion ratio 
100, equilibrium flow). 


Hydrazinium Nitroformate 


HNF is an orange-yellow solid of formula 
N2HsC(NO>)3. HNF is made by a precipitation 
reaction between hydrazine (N2H,) and nitro- 
form [HC(NO2)s3]: 


N>Haq + HC(NO>)30) ~NoHs"*C(NO2)3)_ + 
84 kJ/mol 


As this reaction is exothermic, the process 
vessel must be thoroughly cooled. After reac- 


tion, the HNF crystals are recrystallized, to im- 
prove the purity of the raw HNF and to obtain 
different particle sizes. HNF is now commer- 
cially produced in batch quantities by Aerospace 
Propulsion Products (APP) in Hoogerheide, 
The Netherlands. At the moment, about 50 kg 
of HNF are produced each year. 


Purity is an important factor for the stability 
of HNF. Friction and impact sensitivity increase 
in the presence of even trace impurities. Also, 
the melting point and thermal stability are very 
sensitive to the purity of the HNF. Because of 
this sensitivity, the melting point of HNF is 
used as a simple method to determine its purity. 


HNF particles tend to crystallize into needle- 
shaped crystals, with large length to diameter 
ratio’s (L/D). Needle-shaped HNF is more sensi- 
tive to impact and friction and is not desirable 
for propellant production. Furthermore, the 
rheology of needle-shape-based propellants pre- 
vents the manufacturing of propellants with 
high solid loadings. Cubic or spherical crystals 
would be ideal for both sensitivity and casting 
properties. By recrystallization, the L/D ratio 
can be reduced to values between 4 and 5 (see 
Figure 2), and the different processes also in- 
fluence the mean particle size. At the moment it 
is possible to manufacture HNF crystals with 
particle size varying between 5 and 2000 um. 


Figure 2. Recrystallized HNF with length to diameter ratio’s between 4 and 5, mean diameter about 


500 uum. 
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Monopropellant Combustion 


Because of the highly exothermic HNF de- 
composition, self deflagration of pressed neat 
HNF is possible, at least in the experimental 
pressure range of 0.025 to 10 MPa.'**! The re- 
gression rate, 7», of HNF monopropellants was 
determined by McHale and von Elbe.''®! They 
determined the burning rate of HNF-filled Py- 
rex tubes at approximately 75% of the theoreti- 
cal maximum density (TMD). For safety reasons, 
the HNF was only loosely packed in the Pyrex 
tubes and not pressed. More recently, regression 
rates of pressed HNF at 96% TMD were meas- 
ured.”'"! The results of both measurements are 
shown in Figure 3. Neat HNF regression rates 
are high, reaching more than 40 mm/s at 10 
MPa. HNF also burns nicely at subatmospheric 
pressures down to 0.025 MPa. 
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Figure 3. Experimentally determined 
regression rates of HNF monopropellants 
for two different TMD values. The beginning 
and ending points of the lines indicate the 
experimental pressure range. 


Parr and Hanson-Parr measured the species 
profiles of HNF flames using several optical 
techniques.'''! The transient combustion of HNF 
was studied by Finlinson using a laser recoil 
method.!'7! The HNF oscillations were found to 
be of the same order of magnitude as HMX. 


Decomposition 


To study the decomposition behavior of HNF 
under combustion-like conditions, temperature- 
jump (T-jump) experiments were conducted by 
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Williams and Brill.''’! Species were detected 
with a Fourier Transform Infra-Red (FT-IR) 
apparatus. 


The T-jump studies show a strong depend- 
ence of the species formed during decomposi- 
tion on temperature. It was found that the same 
three temperature regions as for HNF mono- 
propellant combustion can be distinguished: 


1) Decomposition below 123 °C: 
Below the melting temperature of HNF, pre- 
heating only; no decomposition occurs. 


2) Decomposition from 123 to 260 °C: 

This temperature range corresponds with the 
foam zone found during HNF monopropel- 
lant combustion. Ammonium nitroformate 
(ANF, NH,C(NO,)3) aerosol, HC(NO );, 
N>Hy, N2O, H2O, and CO were detected in 
the gaseous phase. HNF aerosol was not de- 
tected; so the evaporation of HNF followed 
by gas decomposition was concluded to be 
negligible. The following reaction scheme 
was proposed: 


2 HNF<) > 2ANF +N, + H+ 
126 kJ/mol 


HNF,(y => HC(NO3)3(¢ + NH, 
— 172 kJ/mol 


2 ANF > NO + 2CO + 4H,0 + 
5/2 O2 + 3 Nz + 360 kJ/mol 


Taken together, the reactions above yield 
314 kJ/mol of HNF. This exothermicity is 
confirmed by the thermocouple experiments; 
this might explain the self-sustained combus- 
tion of HNF at subatmospheric pressures."*! 


3. Decomposition above 260 °C: 

This temperature region corresponds with the 
combustion zone found during HNF mono- 
propellant combustion. Above 260°C, the 
formation of CO, is observed for the first 
time. The amount of CO, increases with in- 
creasing temperature; hence, the exother- 
micity increases with temperature. Above 
350 °C, the amount of ANF, N,O, and CO 
decreases. At 400 °C, the reaction can be de- 
scribed by 


HNF > 2NO + CO; + 2H,0 + 
3/2N. + 1/2 Hy + 623 kJ/mol 
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This reaction is the most strongly exother- 
mic of the above reactions. The products of 
this reaction are not in the thermodynamic 
equilibrium concentrations that exist at the 
theoretical, adiabatic flame temperature. 
Further reaction at higher temperature is 
expected. 


The above reaction pathways differ from 
those proposed by von Elbe et al.'*! By simple 
analytical techniques, they determined that NO, 
was released during HNF monopropellant com- 
bustion. This NO, reacts with the hydrazine 
component of HNF, liberating nitroform, ac- 
cording to 


2 NO, + NH HC(NO3)3 > 
NO + NH4NO; + HC(NO3)3 


Because the free nitroform is unstable well be- 
low 120 °C, this reaction is immediately fol- 
lowed by 


HC(NO,); > x NO, + intermediate products 


where x is presumably >1, but cannot exceed 3. 
It was concluded that at elevated temperatures, 
confinement of HNF in a closed container con- 
stitutes a severe explosion hazard, because a 
run-away increase of the NO, concentration, 
terminating in an explosion, may be expected. 
Effective venting of the container to remove trace 
amounts of NO, would mitigate this hazard. 


The difference between the results from Wil- 
liams and von Elbe is remarkable. During com- 
bustion von Elbe found that NO, plays an im- 
portant role, whereas Williams found no evi- 
dence of NO, in T-jump experiments at com- 
bustion-like temperatures and temperature gra- 
dients. These differences may be explained by 
the fact that NO> is very unstable and might not 
be detected above the platinum filament, be- 
cause it is already reduced to NO by, for exam- 
ple, the released hydrazine. The decomposition 
to NO in one step is very unlikely, for example, 
HMxX and RDX also decompose into NO; first. 


Page 40 


HNF Propellants 


Compatibility of HNF with Propellant 
Ingredients 


Application of HNF in a composite rocket 
propellant requires the use of energetic binders 
for two reasons. As stated earlier, HNF is in- 
compatible with the double bonds in conven- 
tional solid propellant binders like HTPB. Mod- 
em energetic binders do not show this incom- 
patibility as they do not contain double bonds. 
However, apart from the binder prepolymer, 
HNF has also shown incompatibility with many 
isocyanates, which are used for curing the 
binder."'*! A reaction between the hydrazine 
moiety of HNF and the isocyanate groups seems 
most logical, as similar effects were also found 
with hydrazine diperchlorate. 


As HNF is not a very strong oxidizer like AP 
(oxygen balance +21.8 and +34%, respectively) 
performance with hydrocarbon binders like 
HTPB is not maximized due to lack of oxygen 
for binder combustion. Energetic binders have 
energetic groups attached to their backbones 
and often contain oxygen. This reduces the 
amount of oxygen needed to fully oxidize the 
binder decomposition products. Glycidyl azide 
polymer (GAP) is the only energetic binder that 
is produced in large quantities. Therefore GAP 
was chosen as the binder for propellant testing, 
although better performance may be obtained 
with some of the new binders like polyNIMMO 
and polyGLYN, produced by ICI Explosives. 


The azido group in GAP contributes to its 
energy content (Figure 4). During combustion, 
the exothermic scission of the —N3 bond struc- 
ture forms Nx), releasing 685 kJ/mol. Due to 
this highly exothermic decomposition, neat 
GAP gumstock is able to sustain combustion 
with appreciable regression rates (more than 10 
mm/s at 10 MPa).!"*! 


Propellant Combustion 


At the moment there is little experimental 
regression rate data on HNF propellants. TNO- 
PML has tested several propellant combina- 
tions. The most important are summarized in 
Table 2.''*! Due to the crystal-shaped oxidizer 
particles, the high mixture viscosity made casting 
impossible; therefore, the propellants were 
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Figure 4. Molecular structure of glycidial 
azide polymer (GAP). 


pressed. Figure 5 shows the regression rates for 
two different propellants, as measured in a 
chimney burner. The uncatalyzed propellant has 
a very high combustion exponent of 0.81. This 
makes the propellant unsuitable for application 
in rocket engines. By introducing a combustion 
modifier (catalyst), the regression rates are in- 
creased in pressure from 0.1 to 10 MPa. No 
data above 10 MPa are available, but it is ex- 
pected that the catalytic effect will disappear 
above this pressure, and regression rates will 
become similar. Due to the higher regression 
rates at lower pressure, the pressure exponent 
decreased from 0.81 to 0.59. 


Experiments have been conducted to com- 
pare the specific impulse of the HNF/Al/GAP 
propellants to that of conventional AP/Al/HTPB 
propellants.'"*) These experiments have con- 
firmed the theoretical calculations. 


Conclusions 


HNF-based composite rocket propellants 
have a higher performance than any other solid 
rocket propellant at the moment. This, and the 
fact that its exhaust gases are chlorine-free, 
makes it a very likely candidate for replacement 
of AP in many commercial and military solid- 
propellant rocket engines. HNF has a high sensi- 


Modified —~m 


0.1 1 10 
Pressure (MPa) 


Figure 5. Experimentally determined 
regression rates of HNF-based propellants. 


tivity and has shown incompatibility with many 
usual propellant ingredients. Therefore, much 
research still has to be conducted. This makes 
the application of HNF in operative systems 
very unlikely in the short term. However, 
HNF’s advantages are many, and its use as a 
next-generation, solid-propellant oxidizer is 
promising. 
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A Thermodynamic Properties Estimation Method 
for lsocyanates 
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ABSTRACT 


The accurate estimation of the adiabatic 
flame temperature is necessary when evaluating 
a composite propellant. This in turn requires 
additional formulation and properties data in- 
cluding the reactant percentages, chemical for- 
mulas, and formation enthalpies, where the lat- 
ter is sometimes difficult to obtain. When ex- 
perimental data is not available, one must use 
estimation methods. The estimation method or 
its software implementation may not completely 
describe all components or properties of inter- 
est, however. This is often the case with isocy- 
anates. The author has assembled thermody- 
namic data for a variety of isocyanates and, 
from this, has derived Benson Group values that 
may be used in estimating isocyanate properties 
when experimental data is unavailable. Appli- 
cation of these estimates and their effect on the 
estimation of the adiabatic flame temperature is 
considered. 


Keywords: isocyanate, diisocyanate, urethane, 
enthalpy of formation, Benson group, 
estimation, adiabatic flame temperature 


Introduction 


There are many properties of a composite 
propellant that are of interest to the energetics 
chemist. These include, but are not limited to, 
the chemical formulation, adiabatic flame tem- 
perature, and the resulting products gas compo- 
sition and properties as a function of pressure. 
Generally, the chemical composition of the 
propellant ingredients is well established. What 
is sometimes unknown is the formation en- 
thalpy (A;H°) of a particular ingredient. This 
property represents the energy required to form 
this molecule from its constituent elements in 
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their standard states.'! When the combustion 
reaction takes place, and one assumes an adia- 
batic condition, where no heat is gained or lost 
in the system,'*! the enthalpy of the products 
equals that of the reactants. This is the adia- 
batic, or maximum theoretical, flame tempera- 
ture. Thus, knowledge of the reactant formation 
enthalpies is required to accurately estimate the 
flame temperature and combustion properties at 
this temperature. 


Many composite propellants utilize polymers 
to bind together the oxidizer and metal fuel par- 
ticles, and these polymers contribute to the 
combustion process as a fuel source. Polyure- 
thanes are frequently used in this capacity due 
to their relatively low cost, high performance, 
and rheological properties. A polyurethane, put 
quite simply, is a thermoplastic polymer pro- 
duced by the condensation reaction of an isocy- 
anate and a hydroxyl-containing material.”"! An 
isocyanate (NCO) is a nitrogen, carbon, and 
oxygen double bonded together, and this func- 
tional group is bonded to a “base” molecule.'*"! 
A hydroxyl (OH) is a hydrogen bonded to an 
oxygen, which is in turn bonded to a “base” 
molecule.'*"! See Figure 1. 


Acommon propellant binder composition in- 
cludes the long chain hydrocarbon HTPB (hy- 
droxyl-terminated polybutadiene) and the diiso- 
cyanate IPDI (isophorone diisocyanate). The 
diisocyanate differs from the isocyanate in that 
each molecule contains two of the isocyanate 
functional groups. Examples of both iso- and 
diisocyanates are in Figure 2 (located after ref- 
erences). HTPB (see Figure 3) has a hydroxyl 
(OH) group at both ends of the molecule, and 
IPDI has two isocyanate groups (NCO). In the 
process of reaction, the HTPB and IPDI form a 
urethane bond (see Figure 1) that unites them, 
and a synthetic rubber results. 
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Isocyanate plus hydroxyl 


Figure 1. Urethane bond between an 
isocyanate and a hydroxyl. 


Data 


For many isocyanates there is little thermo- 
dynamic data available despite the use of these 
chemicals in tremendous quantities in the plas- 
tics and coatings industries. The urethane con- 
densation reaction mechanism is poorly under- 
stood,”! and the availability of critically evalu- 
ated experimental data is in short supply."! 


A search for thermodynamic data can some- 
times be a lesson in frustration. What is avail- 
able at a conventional library is minimal at best, 
difficult to locate, and often disappointing in 
scope. Fortunately, there are sources of data 


Cis-1,4 20% 


HTPB 


Vinyl-1,2 20% 


that are readily available—ones that provided 
most of the data presented here. 


In particular is an Internet web site main- 
tained by the National Institute of Standards 
and Technology (NIST), NIST Chemistry Web- 
Book."! This online database includes thermo- 
dynamic data for many molecules, of which a 
few are iso- and diisocyanates. As of this writ- 
ing, the ability to search for specific compounds 
is limited to full chemical names, Chemical Ab- 
stract Registry (CAS) numbers, and partial for- 
mulas. 


Another web site, CambridgeSoft’s Chem- 
Finder," has a large database of organic mole- 
cules, and this search engine allows partial 
name searches. When “isocyanate” is used as a 
partial name search, up to 25 “hits” are dis- 
played. Selecting any of these items displays the 
chemical structure (if known), formula, CAS 
registry number, and other properties informa- 
tion. This proved useful in determining the CAS 
numbers for particular isocyanates. The NIST 
site was then used to cross reference the ther- 
modynamic properties when they were avail- 
able. In addition, the structural drawing found 
on ChemFinder proves valuable too, as de- 
scribed later. The 25 hits limit of the Chem- 
Finder web site was eliminated by use of Chem- 
Office Ultra, which includes the complete da- 
tabase found on the web site. 


The data and sources are listed in Appendix 1. 
All data presented is referenced at 298.15K. 
Temperatures are in Kelvin. Units are kJ/mol for 


A;H? (enthalpy), and J/(mol-K) for C,° (heat ca- 


H,cC— CH—— H2C 


n=50 


—CH2 


Trans-1,4 60% 


Figure 3. HTPB (hydroxyl terminated polybutadiene). 
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pacity) and S° (entropy). In the gas phase, only 
A;H° was found; some of these values are esti- 
mates as noted in the data. In the liquid and 
solid phases, only AsH° and C,° were found 
with three exceptions for S°. Some of the data 
was contradictory or apparently erroneous and 
is so noted. Of primary importance in estimating 
the adiabatic flame temperature is the liquid 
phase AH°; the isocyanates and hydroxyl- 
terminated polymers used in propellants are liq- 
uids and react to form a viscoelastic fluid as op- 
posed to a crystalline solid. The majority of the 
A,H?° data found is for the solid phase. 


The Benson Group Additivity Method 


Where experimental data does not exist, one 
must estimate the required properties. There are 
a variety of methods available for the gas phase. 
The interested reader is directed to reference 8 
for lucid descriptions and sample calculations of 
five methods, including that of Benson. For 
composite solid propellant evaluation, the liq- 
uid phase is generally of most interest. The 
method of Benson has been extended to the liq- 
uid and solid phases by Domalski and Hear- 
ing,”! and is available from NIST as the soft- 
ware database program THERM:'""! This pro- 
gram was used extensively in evaluating the 
thermodynamic data presented herein. 


The Benson Group method assigns thermo- 
dynamic values to specific groupings of atoms 
that are commonly found in organic molecules 
containing the elements C, H, N, O, S, Cl, Br, I, 
and F, for which a large body of experimental 
thermodynamic data is available. Specifically, 
the molecule is broken down into groups where 
each group is generally the non-hydrogen atom 
of interest and the atoms immediately bonded to 
it. For example, methane, CHy, is a carbon atom 
with 4 hydrogens bonded to it. In Benson Group 
notation, this is equivalent to C-(H)4. Methanol, 
CH;OH, has 2 groups: a carbon bonded to 3 hy- 
drogens and an oxygen, C—(H);(O), and an 
oxygen bonded to a carbon and a hydrogen, O- 
(H)(C). Carbons that have single, double, triple, 
and resonant bonds are given the notations C, 
Cp, Cr, and Cg, respectively. Unsubstituted 
benzene, for example, would be six C,— 
(H)(Cg)2 groups. There are corrections for a 
variety of molecular configurations. Hexane, 
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CsHi2, when formed into the ring structure 
cyclohexane, imparts additional strain on the 
carbon bonds. The notation used is still six C— 
(H)(C). groups, but then a “correction” group 
“Cyclohexane, RSC” is included to account for 
this strain. When substitutions are made on a 
“base” structure, their proximity to other groups 
can impart additional strain. Examples of this 
include methyl substitution (a “—CH3 correc- 
tion”), and the neighboring interaction ortho 
and meta substitutions, such as o-xylene and m- 
xylene.!'™! 


After all of the component groups are identi- 
fied, their respective contributions are arith- 
metically summed, corrections are applied 
when required, and the resulting values are the 
estimates of the thermodynamic properties. This 
method and its extension to liquid and solid 
phases is surprisingly accurate despite its ap- 
parent simplicity.”! The primary problem in 
evaluating the properties of isocyanates using 
this method is that no values for the NCO group 
have yet been defined"! due to the lack of criti- 
cally evaluated data."*! 


NCO Group Estimation 


The NCO group contributions are required 
to complete the estimation of the isocyanate 
properties. Based on the chemical structures for 
which thermodynamic data was available (see 
Figure 2) only three variations of the NCO 
group are considered based on the specific type 
of carbon-nitrogen bond: a carbon, C, a reso- 
nant carbon Cz, and a doubly bonded carbon 
Cp. The group notations used for these three 
cases are N-(C)(pCO), N-(Cg)(pCO), and N— 
(Cp)(pCO). The introduction of the (pCO) no- 
menclature is to suggest that the CO group is 
itself double bonded to the nitrogen via the car- 
bon rather than using the notation (CO) which 
indicates that two separate CO groups bond to 
the nitrogen. This was also adopted so as not to 
create confusion with the Na and N; nomencla- 
ture used for azo and imino nitrogen bonding. 


Estimating the properties of the NCO group 
is straightforward: using a structural diagram of 
the molecule, like that found on the Chem- 
Finder database, sum the respective contribu- 
tions of all of the other (known) groups in the 
molecule, and then subtract this sum from the 
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cited value. The contribution of the NCO group 
is the remainder. In the case of the molecule 
being a diisocyanate, the remainder is divided 
by two, half for each contributing NCO group. 
Diisocyanates such as IPDI were not used in 
estimating the NCO group values due to the 
molecule having multiple C-NCO bond types. 
The tabulated results are presented in Table 1 
(located after references) and are organized by 
phase—gas, liquid, and solid, and by property 
type—formation enthalpy (A;H°), heat capacity 
(C,°), and entropy (S°). Arithmetic averages 
appear for the various properties and bond 
types. Due to the limited number of data points 
available, no statistical methods were applied to 
the data. One value, as noted, was excluded due 
to its inconsistent and possibly inaccurate value. 
Appendix 2 outlines a brief study of para, meta, 
and ortho group substitutions effects. 


Urethane Reaction Estimate 


One may wish to consider the reaction be- 
tween the component isocyanate and hydroxyl 
compounds that result in a urethane. There are 
no values for the CO-(O)(N) group, how- 
ever.’”'"! The gas and liquid phase A;H° contri- 
butions for this group were estimated from the 
two urethanes listed in Table 2. Due to this very 
sparse data, and for the sake of completeness, 
the author used a different approach to “fill in 
the gap” for the solid phase A;H?, as detailed in 
Table 2. The differences between these single- 
atom estimates and those estimated from the 
two urethanes were 0.82 kJ/mol,,; and —9.33 
kJ/mohjquia for AfH°. This agreement lends con- 
fidence to the solid phase A;H® estimate. These 
are all tentative figures, nonetheless, and should 
be treated as such. 


As an example of this reaction, when a mole 
of HTPB and a mole of IPDI completely react, 


the resulting urethane AsH° does not equal that 
of the reactants due to the rearrangement of the 
atoms taking place in the urethane bond. For 
this particular case, and using the A;H®° values 
in Table 3 combined with the estimated CO- 
(O)(N) group value and bond rearrangements, 
the resulting rubber has a AH° of —230.63 
kJ/mohjiguia 28 Compared to the combined reac- 
tants respective A;H° of 17.65 kJ/mohjiguia—a 
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heat of reaction!”*! of —248.28 kJ/mohjiguia, This 
difference represents two moles of urethane 
bonds due to the reactants having a functional- 
ity of two. Thus, the estimated, exothermic heat 
of reaction per mole of urethane bonds, in this 
case, is —124.14 kJ/mohiiguia. 


Application of Results 


From these tabulated values, one can esti- 
mate the properties for other isocyanates. As an 
example, estimates for IPDI agree reasonably 
well with reference 16 rather than reference 19 
or 21. Appendix 3 compares a sample of esti- 
mated isocyanate enthalpies to values used in 
the propulsion industry. There are significant 
differences in some cases. This leads one to 
consider how the isocyanate AH° estimate im- 
pacts the prediction of the adiabatic flame tem- 
perature of a composite propellant. 


To evaluate this, simulations were run using 
the author’s free energy minimization pro- 
gram.'''! Each simulation was run six times: 
once each for the two cited values for IPDI, and 
once each for the as-estimated A;H° value, this 
value divided by two, this value multiplied by 
two, and for the urethane system described 
above that represents the cured (reacted) isocy- 
anate and hydroxyl-terminated polymer. 


Two generic propellants were used in this 
evaluation, and their formulations and test re- 
sults are presented in Table 3. The diisocyanate 
contributes only a small percentage to the mass 
of the propellant, and thus its overall contribu- 
tion to the propellant enthalpy is minor. Even 
so, the intentionally poor A;H° estimates (times 
2, divided by 2) indicate that there is a subtle, 
but definite, effect on the flame temperature. 
Although it is a small variation, it should not be 
ignored. Accounting for the urethane reaction 
enthalpy resulted in a change to the estimated 
flame temperature that is almost as significant as 
that of the intentional errors in the isocyanate 
enthalpy for these particular formulations. 
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Conclusion 


Experimental thermodynamic data for iso- 
cyanates may be hard to find, but with the data 
that is available, reasonable NCO group values 
may be estimated. These Benson Group NCO 
values accommodate the estimation of some of 
the thermodynamic properties of the arbitrary 
isocyanate, which may then be considered in 
evaluating the composite propellant as a whole. 
The overall effect of the isocyanate enthalpy on 
the adiabatic flame temperature is minimal; even 
so, every effort should be made to acquire or 
estimate accurate AsH° values, especially in 
situations where the effective enthalpy contri- 
bution is greater. The effect of the urethane re- 
action should not be overlooked or ignored ei- 
ther, as it is relatively easy to estimate and ac- 
count for using similar group additivity meth- 
ods. 
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Figure 2. Isocyanate and urethane molecular structures. 
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10 Toluene diisocyanate 11 4,4’-Biphenyldiisocyanate 12 Benzene, 
26471-62-5 2761-22-0 1-chloro-3-isocyanato- 
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Figure 2. Continued. 
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Table 1. NCO Benson Group Contribution Estimates. 


N-(Cg)(pCO) 


Fig. 2 a 
ee. No. CAS # Coo 


-71.98 | 63.94 | NIA | -77.76 | 55.31 


Units A;H° kJ/mol 
C,° kJ/(mol*K) 
s° kJ/(mol*K) 


N-(C)(pCO) 


= 109-0050 1) 70.04 | 80 CT os coca 
114-364 


5124—30-1 i 
Average: : -42.12 | 55.74 } 112.91 | -72.74 | 45.99 


N-(Cp)(pCO) 


3555-94-0 | -40.32 


[a] Value is possibly not accurate and has been excluded from the average. 


[b] These gas phase enthalpies are based on estimates. 
[c] CAS 28182-81-2, 1,6-Hexamethylene diisocyanate CgH).N2Oo. 


Page 50 Journal of Pyrotechnics, Issue No. 6, Winter, 1997 


Table 2. CO-(O)(N) — Group Estimates. 


Estimated Values for the Group CO-(O)(N) Experimental Data 


The resulting CO—(O)(N) contribution: 


CAS # Gas (AH?) Liquid (AH?) Solid (C,°) 


Average: 


Known Benson Group/Single-Atom Contribution Estimates 


The following THERM Groups were used to derive the average contribution of a single-atom to a CO 
(carboxyl) group for A;H°: 


CAS # Gas (A;H°) Liquid (A;H°) Solid (AH?) 


Group 116, CO-(O), 
Group 289, CO-(N)z2 


Dividing each of the above by 2 gives the single-atom contribution. Adding these single-atom values, 
then subtracting this from the known value for the actual group gives the methods residual error: 


—137.24 —149.37 —153.60 
—122.28 -137.38 —140.48 
— 14.96 — 11.99 — 13.12 
—133.26 —185.00 —194.60 
—121.34 -171.63 —180.53 
— 11.92 — 13.37 — 14.07 


The error is evenly distributed, and for all phases it is —13.24. 


The atoms are added together, and the average error added to that sum to yield values for an estimated 
Group CO-(O)(N) A;H?°: 


CO-(O)(N) Group A;H?°: 
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Table 3. Isocyanate A;H° Estimate Effect on Adiabatic Flame Temperature. 
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Application of NCO Group Estimates 


Propellant Ingredients 


ingredient 


HTPB Ca99H39202 


2738.60 
IPDI C12HigN202 222.29 
Ammonium Perchlorate NH,CIO, 117.49 
Aluminum Al 26.98 


Propellant Compositions 


365.22 kJ/mol!*! 
variable 
-295.3"! 

0 


Propellant One Propellant Two 


Estimate / 2 


Value from reference 18b. 


In both propellants, the NCO:OH ratio is 1:1. 


Value estimated by the author using NIST THERM (reference 10), based on structure in reference 22. See 


Flame temperature is estimated using reference 11 and convergence is computed using a binary subdivision 
between 298.15 K and 6000 K, until the difference between successive iterations is less than 1.0 degree K 
as the reactant minus product enthalpies approaches zero. 


The urethane is equivalent to the estimated HTPB and IPDI reacting to form a urethane, and the resulting 
product estimate includes the CO-(O)(N) estimated A;H°. The value used in the simulation for the resulting 
product has the formula C;7H329N2O4, molecular weight 2960.88, and A;H° of —230.63 kJ/mol. 
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Appendix 1 


86—84—0 | Naphthalene, 1-isocyanato- 


“Solid “AH? = 
101-68-8 | Benzene, 1,1'-methylenebis [4-isocyanato- 
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102-36-3 | Benzene, 1,2-dichloro-4-isocyanato- 
ares eorcs ge a armen et ee 


103-—71-9 | Benzene, isocyanato- 


104—12-1 | Benzene, 1-chloro-4-isocyanato- 

sears CON oneal el cells (Oia el a tueatortnamncunaed 
104—49-4 | Benzene, 1,4-diisocyanato- 

See ee 1 ge ge ee 


109-90-0 
Sachaainel festa eed 


eee ey 


112-96-9 
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‘Sold ics Re SARA SGT 
3555-94-0 | 

Siciadt tnnceanetnamaeel joe ee 
4098-71 ~9 | 


51 24-30-1 
ee eal LOIN 2) li SNS MN eaten cantata 
~ 440.6 


‘[** Phase unknown 
26471-62-5 T -84- 


*“ Phase unknown 
96775-58—7 2- “Methyl- el 5- -naphthalene diisocyanate 

cacehe aaa cca a a nl ala 
64711-83-7 | Ethyl-m-phenyle 


2761-22 0 | 
2909-38 8 | 
3173-72 6 | 
3320-—83-0 | Benzene, 1-chloro-2-isocyanato- 


[Mixture sd 
68239-06-5 | Di i 


= [1236.33 0 | 16 | ** Phase unknown 
—872.60 ** Phase unknown 
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Appendix 2 


The groups used in estimating the NCO contributions included corrections where required, such as 
ortho CI-Cl, meta, CH; tertiary and quaternary substitutions, etc. The available experimental data is 
too limited to draw any conclusions in terms of the development of additional correction groups for 
NCO-NCO, CH;—NCO, and CI-NCO next-to-nearest neighbor interactions. Nonetheless, MOPAC [12] 
was used to investigate these interactions. The method used was to draw the structure in ChemDraw 
Pro,"?! import it into Chem 3D Pro,''*! and then minimize the strain of the structure using the PM3 
Hamiltonian. 


The two structures studied were benzene and the chair conformation of cyclohexane. An NCO 
group was substituted, and a second group X was substituted in the para, meta, or ortho position rela- 
tive to the first group. This second group X was either CH3, Cl, or NCO. The results are tabulated be- 
low in both A;H®,,; and the difference between the para value and that of the current position. The in- 
teraction did not become really significant until the ortho position, and in most of the available data, 
this condition does not arise. 


The exception to this, for which A;H®,,,; is available, is a meta NCO—NCO interaction, 64711—83-7, 
ethyl-m-phenylenediisocyanate. The cited A;H®,,; is -135 kJ/mole. The MOPAC estimate is —47.18. 
This represents a significant difference that cannot be easily evaluated without additional gas phase 
enthalpy data. 


Benzene 


ete HM 8 Ne 


Cyclohexane 


a 
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Appendix 3 


Comparison of Estimated and Cited Isocyanate Enthalpies 


Fig. 2 Estimated A;H° Cited 
26 


822-06-0 1,6-Hexamethylene diisocyanate 


—289.26 —248.76 —331.12 —226.18 
—504.70 


4098-71-9 Isophorone diisocyanate 


—335.15 —347.57 —398.33 —371.62 
—466.21 
—466.07 


15646-96-5 2,4,4 Trimethylhexamethylene diisocyanate 
—362.15 —332.22 —414.13 —377.61 
26471-62-5, Toluene diisocyanate 


—104.16 —141.67 —142.27 —24.05 
—623.17 


68239-06-5 Dimeryl diisocyanate 


—871.45 —972.86 —1125.06 —1236.33 
—872.60 


1) The phase of the cited values listed is unknown. 


2) The chemical name of the cited value, and its value, were all that was available for the above. It is assumed 
to be the same as those presented by CAS registry number. 


Appendix 4 


Suggested Values for the NIST THERM Database file THERM.COD 


The Group Numbers chosen, 800-805, were an arbitrary choice by the author. Please note that the 
user takes full responsibility in editing this file. 


800!N-(CB)(DCO)!-62.50!*!*!-71.98163.94!*!-77.76155.31!106.90!! 
801!N-(C)(DCO)!-75.07!*!*!-42.12155.7411 12.91!-72.74145.99151.97!! 


805!CO-(N)(O)!-125.22!*1"1-159.571*!*!-176.70!188.56!"1! 
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A Severe Human ESD Model for Safety and 
High Reliability System Qualification Testing 


Richard J. Fisher 
Formerly with Sandia National Laboratories, Albuquerque, NM 87185, USA 


ABSTRACT 


A severe human electrostatic discharge 
(ESD) equivalent source model has been devel- 
oped for use in qualification testing of systems 
that have stringent safety or very high reliabil- 
ity requirements. The model produces the now- 
acknowledged worst-case waveform, and the 
values of the constituent components have been 
selected from measured human body electrical 
parameters to maximize the peak amplitude and 
rate of rise of the short-circuit discharge cur- 
rent and energy transfer to the victim system. 


Introduction 


Many systems have either stringent safety or 
very high functional reliability requirements. 
Often these same systems, military ones in par- 
ticular, are subject to operations and mainte- 
nance in field or repair depot environments 
where there is no guarantee of ESD source con- 
trol. This combination of requirements implies 
the need for a capability to demonstrate system 
tolerance to severe human ESD. 


Most of the broadly accepted existing hu- 
man ESD simulation circuits, for example those 
in the MIL-STD-883C, International Electro- 
technical Commission (IEC), and European 
Computer Manufacturers Association (ECMA) 
models, appear to incorporate average values 
for equivalent body electrical parameters and 
initial charge voltages. They therefore corre- 
spondingly replicate typical, or moderate, ESD 
events. While these environments adequately 
address the testing needs for a large variety of 
commercial and other systems, they are inade- 
quate for qualification of systems of the type 
mentioned above. 
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In the present effort, a database of the better 
documented measurements of human body elec- 
trical parameters—capacitance, resistance, and 
inductance—and the distribution of acquirable 
body voltages was assembled from the litera- 
ture. The data were found to be surprisingly 
sparse, especially with respect to body voltage. 
A severe human ESD source model was devel- 
oped with component values selected from the 
available data to yield short-circuit discharge 
parameters with reasonable upper-bound values 
that could be encountered by a system over its 
life-cycle under actual conditions. 


Elimination of the attachment arc in ESD test- 
ing is now widely recognized to enhance both 
the repeatability of test results and their correla- 
tion to failure distributions observed on the 
same test objects due to actual human ESD.!"7! 
Therefore, the present model is intended for di- 
rect connection to the test item, and no elements 
corresponding to an arc are included. 


In this paper, the rationale is given for the 
form of the model and its component values. For 
brevity hereafter, the term ESD is used to mean 
human body ESD unless otherwise indicated. 


Equivalent Source Model 


Output Waveform 


The key aspects of pulsed electrical over- 
stress transients in general, and ESD in particu- 
lar, that influence the nature and level of in- 
duced effects are peak amplitude, rise rate, and 
total energy transferred to the victim object. 
The correlations between effect and the peak 
amplitude or total energy deposited by the pulse 
can be readily appreciated. The extreme sensi- 
tivity of victim systems (both electronic and 


Page 57 


ordnance) to rise rate of the discharge current is 
perhaps less intuitive, but nonetheless well- 
documented [e.g., Refs. 1-5]. Effects or failure 
mechanisms that are sensitive to the derivative 
of the driving current include radiation from the 
source current, magnetic field coupling into 
internal wiring or circuit boards, inductive volt- 
age transients, and others. 


ESD events result in a very wide variety of 
waveshapes. It is now widely accepted, how- 
ever, that the worst-case waveform is of the 
type shown in Figure 1, which corresponds to 
the short-circuit current often observed in dis- 
charges from a hand grasping some small me- 
tallic object. That this particular shape repre- 
sents the worst case follows from the presence 
of the initial narrow spike, which generically 
reflects the highest peak amplitudes and fastest 
rise rates observed in the extensive empirical 
database reported by King, Richman, and others 
in recent years.*”! This waveform was there- 
fore adopted as the desired output of the model. 


4ns/div 


Figure 1. Worst-case human ESD waveform 
(After Richman!”). 


Source Circuit Model 


The traditional human body ESD equivalent 
source models have been simple series RLC 
circuits in which the R corresponds to the resis- 
tance of that portion of the body in series with 
the discharge, L is the similar body inductance, 
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and C is the capacitance of the body to its sur- 
roundings. The body inductance is often ne- 
glected, as, for example, in MIL-STD-883C, in 
apparent tacit recognition of the fact that zero 
inductance represents the worst case. 


Recent refinements to source models have 
been suggested to account for various observed 
features in actual ESD currents, such as multi- 
ple-pulse discharges"! and the important initial 
narrow spike of Figure 1°! Multiple pulses oc- 
cur only in discharges involving sparks. Since 
total energy transferred is the sum of that avail- 
able to the victim in the individual pulses, 
worst-case stress can be simulated in a single 
current pulse from an appropriate source. Com- 
parison of the minimum interval between indi- 
vidual pulses (~20 ps'*!) versus the thermal time 
constants of victim electronics (hundreds of 
nanoseconds or more) also indicates that the 
transfer of all available energy from the source 
in a single pulse represents the worst case. 


All ESD source models rely on the Thevenin 
linear circuit equivalency principle, which 
guarantees that if the open circuit voltage and 
short-circuit currents are identical at the output 
of any two circuits, then the two circuits are 
completely equivalent electrical sources. The 
validity of test results obtained with this type or 
model therefore corresponds in part to the de- 
gree to which the ESD event being simulated is 
linear. In that regard, there are two potential 
sources of nonlinearity with respect to initial 
charge voltage to be considered: equivalent body 
electrical parameters and pre-breakdown corona 
and other effects in spark discharges. No solid 
evidence was found in the literature that there is 
any significant nonlinearity associated with the 
former, although there appeared to be a hint in 
one study of decreasing dynamic body resis- 
tance with increasing initial voltage.''®! Most of 
the available data are too sparse and involve too 
many simultaneous variable changes to allow 
analysis for that particular effect. Complica- 
tions with respect to reduction of peak current 
and rise rate due to corona effects!!" are 
avoided altogether in the present model, since it 
is intended for direct connection to the test item 
prior to discharge. 
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Figure 2. Dual-RLC equivalent human ESD circuit model due to Richman." 


The basic form of Richman’s 1985 model 
(Figure 2), shown by him to be capable of re- 
producing the important early spike,'! was 
adopted as the preferred baseline. Here, the 
components labeled Cy, Ry, and Ly are intended 
to correspond physically to the capacitance and 
other electrical parameters of the arm and hand 
as they reach out towards the victim. From this 
viewpoint, Figure | can be interpreted as a com- 
posite of a very fast and narrow initial pulse of 
limited energy content, corresponding to the 
rapid discharge of the hand capacitance, super- 
imposed on the slower, broader discharge asso- 
ciated with the bulk body capacitance. 


0.1 pH 


The present model (Figure 3) represents a 
slight modification of the form of Richman’s 
circuit, in that the resistance and inductance 
associated with the hand and arm have been 
moved from a parallel branch into series with 
the discharge path to improve correspondence 
with physical reality. In the following, the data- 
base on equivalent body electrical parameters is 
reviewed, and the basis for selection of the spe- 
cific element values reflected in Figure 3 is 
given. 


1100 


Figure 3. Modified human ESD source model for severe ESD simulation. 
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Electrical Parameters of the Human Body 


Although values for body parameters are 
widely asserted in numerous ESD handbooks 
and papers, references in the literature actually 
describing either the data themselves or the 
methodology under which they were obtained 
are surprisingly meager. On the other hand, in 
the few better documented cases that were 
found, there is a rather comforting degree of 
agreement among data that were independently 
acquired using a considerable range of techni- 
ques.""!?"41 Among these were the direct use 
of commercial capacitance bridges and various 
schemes for extracting body R and C from re- 
corded discharges from volunteers charged to 
controlled initial voltages. Capacitance data are 
reported for subjects of differing physical char- 
acteristics (height, weight, shoe size, etc.) and 
for various shoe sole and flooring material 
combinations. Table 1 summarizes these results. 
The extreme upper and lower values noted were 
105 and 1100 pF, respectively, although several 
other authors alluded to body capacitances as 
low as a few tens of picofarads without descrip- 
tion of the source of data. Sullivan and Under- 
wood!'*! provide fitted normal distributions for 
their data (Figure 4). The low-end skirts of these 


Table 1. Mean Measured Body Capacitance. 


Standing Sitting 
Investigator (pF) (pF) 


Sullivan and 
[Sperber and Blink’ | 130 | 300. 
Calvin et al.!'*! 


164/67 


curves indicate values of this order, although 
nothing that low appeared explicitly in their 
tabulated results. Of more practical import, 
however, is the fact that, according to their dis- 
tributions, nothing over 400 to 450 pF would be 
expected, even in situ in an automobile envi- 
ronment, which intuitively might be expected to 
yield the highest capacitance due to the increase 
in body surface area in close proximity to its 
surroundings. 


SITTING 
IN 
AUTOMOBILE 


(pF) 


Figure 4. Measured body capacitance distributions for individuals under in-plant and automobile 


environments (After Sullivan and Underwood'"*'). 
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Body resistance apparently exhibits a much 
broader range of values, as is suggested by the 
data in Table 2. It is most dependent on skin 
moisture, the particular portion of the body in- 
volved in the discharge path, and the contact 
area and degree of pressure applied in grasping 
a metal object when some hand-held article is 
included. A distinction is made in Table 2 be- 
tween the two main types of discharges from 
the hand, namely, those with and without the in- 
volvement of a metal object. Qualitatively, all of 
the data reviewed indicate that the latter case 
implies the highest effective resistance. Other 
values appearing in the literature without dis- 
cussion of their origin range from 100 to 2500 
ohms. The only values for body inductance 
found were between 0.1 and 2 uH. From this 
database, values for the various circuit elements 
were selected according to the following ration- 
ale. 


Table 2. Measured Effective Human Body 
Resistance. 


Fingertip | Metallic Object 
Investigator (ohms) (ohms) 
5 


1713/308| 535/172 


First, in order to maximize the available dis- 
charge energy, the highest reasonable bulk body 
capacitance was selected. While the single high- 
est value reported in any of the better docu- 
mented studies reviewed was 1100 pF,'"*! this 
value was a factor of 3.5 higher than the next 
lower one (314 pF). Furthermore, from Fig- 
ure 4, it is seen that a value of 400 pF represents 
something like the 40 point of the Sullivan and 
Underwood data. Hence, in the spirit of a practi- 
cal extreme, a value of 400 pF was selected for 
the main body capacitance Cg of the present 
model. 
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Total body resistance is the sum of that as- 
sociated with the main trunk and that due to the 
arm and hand. In order to maximize both peak 
current and energy transfer to the victim, source 
resistance was minimized by selecting the low- 
est reasonable value of total body resistance 
found during the review. This is Tucker’s 357 
ohms (Table 2), which was rounded up to 360 
ohms. This total was allocated between the main 
body (Rg) and the arm and hand (Ry), 250 ohms 
to the former and 110 to the latter. 


Hand capacitance Cy was chosen to be the 
upper limit of Richman’s range"! to maximize 
the amplitude and energy in the initial spike. 
Hand inductance Ly was then selected to be 0.1 
uH to provide a 10-90% risetime for the initial 
peak of just under a nanosecond, corresponding 
to the faster risetimes evidenced in the data of 
King and Reynolds and Richman. Finally, a 
value of 0.5 wH for the less critical main body 
inductance Lz was adopted, also from the range 
of values reported by Richman. 


It should be noted that both the severe 
waveshape definition and its equivalent model 
are based on actual short-circuit currents pro- 
duced in configurations with very low return 
path inductance.?*7! Although other discharge 
scenarios could, in fact, involve substantial re- 
turn path inductance, which would significantly 
affect the discharge current, this will not neces- 
sarily be the case. Hence, once again in the 
spirit of conservatism, no extraneous inductance 
is incorporated into the source model itself. 


Initial Charge Voltage 


Very few well-documented data were found 
regarding the ranges of electrostatic voltages 
acquirable by personnel under various envi- 
ronmental circumstances. In principle, the abso- 
lute upper limit on the voltage that can be sus- 
tained is set by the point at which charge leak- 
age mechanisms, primarily corona, provide equi- 
librium with respect to further charge accumula- 
tion. This implies a practical upper limit on a 
person in normal proximity to his surroundings, 
variously stated to be in excess of 25 kV"! and 
30 to 39 kv 4617 
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Table 3. Body Voltages and Capacitances after Exiting from an Automobile. 


{14] 


Voltage Shoe Type 


3/8-inch Plastic Intermediate Size 
3/8-inch Plastic Station Wagon 
3/8-inch Plastic 

3/8-inch Plastic 

3/8-inch Plastic 

3/8-inch Plastic 


Rubbers Over Shoes 
5/8-inch Rubber Sole 
3/16-inch Leather Sole 


3/4-inch Plastic Sole 
3/4-inch Plastic Sole 


Sullivan and Underwood!'"*! report measured 
data from a sample of subjects within labora- 
tory and production plant environments. They 
give mean values of 5.8 kV and 1.4 kV (no 
standard deviations) for low and high humidity 
conditions, respectively, but point out that their 
data are skewed in the low direction due to 
limitations of their measuring equipment. Sulli- 
van has privately reported measured personnel 
voltages in excess of 70 kV within an airplane 
environment.''*! Sperber and Blink''*! give 
measured data (Table 3) on a sample of people 
with varying heights and weights as they exited 
from automobiles. Their data contain the high- 
est measured personnel voltages reported among 


140 
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a | 
© 100 f) 
oO | 
= 80 
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c 60 | 
8) 
= 40 
3S 
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the better documented studies. Conspicuously 
absent from the database appear to be any 
measurements of personnel voltages acquired 
under outdoor conditions and activities. If lev- 
els of 20 kV or more can be developed under 
plant-type condition, one can readily envision 
significantly higher voltages in the case of, say, 
a soldier or airman operating outdoors in a dry 
wind. In view or these considerations, a reason- 
able upper bound for use in the present model 
was chosen to be 25 kV, although arguments 
could be entertained for an even higher value. 


20 0 50 100 150 200 250 300 
TIME (ns) 


(b) 


Figure 5. Output of the ESD model of Figure 3 computed for an initial charge voltage of 25 kV; 


(a) early and (b) longer time base. 
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Model Output 


Figure 5 gives the short-circuit output cur- 
rent computed with the model shown in Figure 3 
for an initial voltage of 25 kV. The waveshape 
compares favorably with the generic worst-case 
waveform represented in Figure 1. Although the 
model was specifically developed for use in 
system-level qualification testing at severe lev- 
els, it contains only linear elements; and its out- 
put amplitude therefore scales linearly with ini- 
tial voltage. This allows it to be used at lower 
voltages while still preserving the desired wave- 
form. The dashed early portion of Figure 5a 
represents the component or current due to the 
main body capacitance. Of the three main stress 
parameters, peak amplitude, rise rate, and total 
pulse energy, the initial peak accounts for the 
first two. The energy contained in the first peak, 
relative to that in the sustained current, is ap- 
proximately the ratio of the capacitance of the 
hand to that of the body, which in this case is 
0.025. 


Conclusion 


The ESD model presented herein provides a 
short-circuit output current corresponding to the 
acknowledged worst-case waveform evidenced 
in the available human ESD empirical database. 
Values of the constituent electrical components 
were chosen from the ranges of substantiated 
measurements of human body equivalent elec- 
trical parameters and acquired charge voltages 
reported in the literature. The severity of the 
resultant output is thought to replicate a reason- 
able worst-case ESD transient in terms of rise 
rate, peak amplitude, and total pulse energy 
available to the victim. The model was devel- 
oped for specific application as a human ESD 
standard at Sandia National Laboratories for 
systems required to tolerate severe ESD. How- 
ever, since its output is linear with initial volt- 
age, the model can be used at more moderate 
levels as well. 
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Computational Mechanics Assoc. Courses: 
Material Behavior at High Strain Rates 
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Detonation Physics 
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Conact: Computational Mechanics Associates 
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Consideration of Alternate 
Whistle Fuels 


Rembert Amons 
Veldheimlaan 23, 3702 TA Zeist, 


The Netherlands 
e-mail: rembie.amons@tip.nl 


ABSTRACT 


The following article aims to be a basic ac- 
count of the use of phthalic acid salts in whis- 
tle compositions. The experiments I conducted 
with these compositions began when I came 
across the comments on modern Chinese whis- 
tle compositions in the new edition of Ronald 
Lancaster’s Fireworks Principals and _Prac- 
tice. The writer reported problems with the use 
of these substances which stimulated me to 
investigate the matter in practice. The first part 
of this article deals with the properties of the 
basic materials, the second part is an experi- 
mental report and points at a few requirements 
for the successful use of these salts in whistle 
mixes. Finally some advantages of these salts 
over the more common materials are discussed. 


Benzoic Acid 
m.p. 121 °C 


Journal of Pyrotechnics, Issue No. 6, Winter, 1997 


Phthalic Acid 
m.p. 210 — 211 °C 


The Basic Materials 


Modern whistle compositions typically con- 
tain a salt from an aromatic carboxylic acid as 
a fuel in combination with potassium chlorate 
or potassium perchlorate as the oxidizer. Salts 
like potassium benzoate and sodium salicylate 
have found extensive use and seem by far the 
most popular choice judged by their prominent 
presence in the literature. However, some refer- 
ence works over the years have mentioned other 
substances which are used in fireworks from 
the East.!"! 


The alternatives mentioned are the salts 
from phthalic acid or benzene dicarboxylic 
acid. Basically the difference from benzoic 
acid, which may be used to produce potassium 
benzoate, is that phthalic acid is a polybasic 
benzene carboxylic acid. That is, it has two 
carboxylic acid groups instead of one. The situ- 
ation with the phthalates is therefore a bit more 
complicated than with the benzoates. 


First of all, phthalic acid or benzene dicar- 
boxylic acid, comes in three different isomers 
which do not necessarily have the same prop- 
erties. In fact in some respects they do have 
different properties, and so do their salts. 


Secondly, phthalic acid forms salts in more 
than one way and either one or both acid 
groups may react with a base. The three differ- 
ent isomers are: 1,2 benzene dicarboxylic acid 


COOH 


COOH 


Terephthalic Acid 
m.p. >300 °C, sublimes 
without melting 
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(the ortho isomer; phthalic acid), 1,3 benzene 
dicarboxylic acid (the meta isomer; isophthalic 
acid) and 1,4 benzene dicarboxylic acid (the 
para isomer; terephthalic acid). The author has 
not experimented with the meta isomer yet, but 
the costs associated with this material might 
render any results of academic interest only. 
This leaves us with phthalic acid and 
terephthalic acid for the moment. Although 
these substances share the same formula, the 
molecule looks quite different. The para posi- 
tion of the carboxylic acid groups typically 
give a very low solubility in water, whereas the 
ortho position gives a good solubility. These 
properties must be reckoned with when pro- 
ducing the salts from the respective acids by 
neutralization with a base. 


It is by no means be the only possibility, but 
the potassium salts seem the most likely candi- 
dates as whistle fuels. One of the salts that has 
been mentioned in the fireworks literature is 
potassium hydrogen phthalate. This is the salt 
whereby one acid group is neutralized. It is 
sometimes referred to as KHP, and it is readily 
soluble in water (107°°/33'°° grams in 100 cc). 
Because of the good solubility of phthalic acid, 
this salt can be comfortably produced with po- 
tassium carbonate. Potassium hydrogen tere- 
phthalate on the contrary is only sparingly 
soluble in hot water. It is said to be used exten- 
sively by the Chinese.”! In this case the 
acid/base neutralization is more difficult be- 
cause of the low solubility of terephthalic acid 
and the reaction product. It can be done, but 
care must be taken that the reaction is com- 
plete. Both salts are suitable as a fuel in whistle 
mixes, but it must be noted that the salt from 
the para isomer is less likely to absorb moisture. 


A salt of different properties will be formed 
by complete neutralization of the acid. In this 
case only the salt from the para isomer, dipotas- 
sium terephthalate, is usable. The salt from the 
ortho isomer is useless because it absorbs 
moisture strongly. But even with the salt of the 
para isomer, it is recommended to allow for a 
slight excess of acid to reduce moisture prob- 
lems. The production of this salt is easier be- 
cause the reaction product is readily soluble. 
This salt makes a good whistle fuel but has 
more affinity towards moisture than the hydro- 
gen phthalates (about the same as sodium sali- 
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cylate). Small amounts of moisture will have a 
significant negative effect on performance.""! 


Whistles with Phthalic Acid Salts 


As stated, both KHP and potassium hydro- 
gen terephthalate will work in whistle compo- 
sitions. They can be used in much the same way 
as potassium benzoate or sodium salicylate, 
though the performance can be somewhat less 
than with the former in terms of whistle power. 
Some, however, seem to be able to make their 
whistle rockets ascend to greater heights with 
the use of KHP."*! 


To make these compositions perform close 
to the level of the benzoate based compositions 
a few important requirements have to be met. 
The first important point is the particle size. 
Normally materials passing #120 British Stan- 
dard (BS) Sieve are recommended for whistle 
mixes. This will not suffice for the hydrogen 
phthalates if the average particle size is not 
way below that. Tubes (8, 10 mm ID) charged 
with potassium perchlorate and potassium hy- 
drogen (tere)phthalate both passing a #120 BS 
sieve burned in a sputtering manner with short 
eruptions of whistle sound. This tendency 
seems to be strongest with whistle mixes with 
an oxidizer/fuel ratio around the stoichiometric 
point (for example, 2.3:1). Compositions with 
an excess of oxidizer (3:1) do not show this 
tendency as much, but the whistle sound seems 
somewhat softer. The irregular burning prob- 
lem can be mitigated by the addition of at least 
1% catalyst such as iron(III) oxide or copper(II) 
oxychloride (note: it is unwise to use this cata- 
lyst with potassium chlorate). It must be said 
that whistle mixes made with these salts bene- 
fit clearly from a catalyst. Further reduction of 
the particle size, if possible below #300 BS, 
also improves performance, but a catalyst might 
still be necessary for an optimum whistle sound. 


A modification of the burning front such as 
a progressive burn taper will also help. In fact, 
superior whistles can be made this way. A 12 
mm ID whistle rocket with about 15 mm of 
free space over the burning front and a cylin- 
drical core of 6 mm ID and 15 mm in composi- 
tion length will take a 20-gram load to an ap- 
propriate height using about 5 grams of whistle 
composition. The composition is a potassium 
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perchlorate/potassium hydrogen terephthalate 
mix in the ratio 2.3:1 with 1% additional cata- 
lyst [iron(II) oxide, Fe.03]. The rocket is bot- 
tom lit (although it will not explode when core 
lit) and ascends with an aggressive how] that is 
quite distinct with this fuel. The dimensions of 
this whistle motor are adjusted to combine a 
high level of thrust and sound. Enlarging the 
burning surface through extended core length 
will increase the thrust, but not necessarily the 
sound effect. The pressure in the burning cham- 
ber of the whistle device should not be too high, 
or it will suppress the sound output. This phe- 
nomenon is known among whistle rocket 
builders, although the mechanism behind this 
seems not entirely clear. 


Whistle Composition Casing 


The second important point is the loading 
pressure. Consolidation with hand pressure 
only, has proved to be unreliable and will re- 
duce the whistle power of these mixes. In fact 
the loading pressure needs to be higher than 
with the mixtures based on potassium benzoate 
to obtain similar results.!*! Frequently the whis- 
tle will burn without whistling right away if 
the composition is consolidated with insuffi- 
cient pressure. 


Possible Advantages of Phthalic Acid Salts 


This is perhaps the most important question 
of this article: do these new materials have any 
advantage over the more common ones? The 
answer is both yes and no. Let's start with the 
drawback. The cost of these materials will 
probably be higher than the cost of the benzo- 
ates or salicylates, although everything de- 
pends on finding the right source (some sug- 
gested the plastics industry and this might in- 
deed be a possibility). 
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The advantages are multiple. First of all, 
the low hygroscopicity of potassium hydrogen 
terephthalate will make this fuel an excellent 
choice in damp conditions. Potassium hydro- 
gen terephthalate is also easy to reduce to a 
fine powder (potassium hydrogen phthalate is 
a bit more difficult in this respect). But more 
importantly, tests have shown the sensitivity to 
friction of the basic potassium perchlorate/po- 
tassium benzoate mix to be twice as high as 
with the potassium perchlorate/potassium hy- 
drogen terephthalate mixture.! Needless to 
say, an additional catalyst might increase the 
sensitivity. 


Of course these matters should be investi- 
gated more extensively and with proper refer- 
ences. The information contained in this article 
is therefore to be regarded as only an indica- 
tion. In practice, however the whistle mixes 
based on the phthalic acid salts perform well 
and have proven to be reliable. Hopefully this 
new material will give rise to interesting ex- 
periments. 
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Frequency Stabilization of 
Large Diameter Strobe 
Flares 


K. L. Kosanke 
PyroLabs, Inc., Whitewater, CO 81527, USA 


Pyrotechnic strobe devices burn to produce 
a series of approximately equal intensity light 
pulses at approximately equal time intervals, as 
illustrated in Figure 1. The typical pyrotechnic 
strobe process"'! begins with a dark reaction 
during which a reactive slag is produced. Subse- 
quently, hot spots develop within the slag, 
eventually triggering a flash reaction which 
consumes most of the slag to produce a bright 
flash of light. The process continues with the 
dark reaction again producing more reactive 
slag, which then eventually flashes. The proc- 
ess continues in this way to produce a series of 
short bright light flashes with relatively long 
periods of darkness in between. This pyrotech- 
nic strobe process is illustrated in Figure 2. 


There is interest using large diameter pyro- 
technic strobe flares as signals and for potential 
military application.'”! However, as the size of a 
strobe flare is increased, there is a tendency for 
the frequency of the flashes to increase and 
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Figure 1. An illustration of the light output 
from a pyrotechnic strobe device. 


become erratic in period and intensity. The 
result tends to more of a random flickering, or 
even continuous light output, than a true strobe 
effect with clearly defined dark periods be- 
tween the flashes. This brief article has been 
prepared to suggest two potential methods to 
help control the burning of large diameter 
strobe devices, for which the author does not 
have time to investigate. 


The idea for stabilizing the burning of large 
pyrotechnic strobes is analogous to a common 
method used to control analog electronic oscil- 
lators. It is often the case that a free running 
oscillator, with a natural frequency a little 
lower than desired, is controlled and stabilized 
by applying regular synchronization pulses, a 
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Figure 2. An illustration of the mechanism of pyrotechnic strobe burning. 
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little before each natural oscillator cycle is 
completed. In a somewhat similar process, it 
may be possible to control large strobe de- 
vices. If a strobe composition is formulated 
that has a low natural flash frequency, upon 
burning it will have prolonged periods of dark 
reaction during which slag is produced. (In 
fact, it might even be preferred to use a com- 
position that produced the needed reactive 
slag, but, which never produced the hot spots 
needed to trigger the flash reaction.) Then, if at 
regular intervals the energy needed to trigger 
the flash reaction is supplied from another 
source, light flashes might be produced with 
the same regular intervals. 


Two methods for providing the periodic en- 
ergy for the flash reaction, have been consid- 
ered. One method might be to incorporate a 
relatively small diameter core of normal strobe 
composition into the large diameter strobe de- 
vice, such as illustrated in Figure 3. In that way, 
the energy regularly produced as the small 
strobe core flashes, might be used to trigger 
controlled flash reactions of the slag across the 
surface of the large device. 


Low Frequency 
Strobe Composition 


Core of Higher 
Frequency Strobe 
Composition 


Figure 3. Illustration of a small core of strobe 
composition used to stabilize a larger strobe 
device. 


In Figure 4, an alternative flash triggering 
energy source is suggested. If the large strobe 
device has a pair of electrical conductors in- 
corporated into it, it might be possible to pro- 
vide the triggering flash energy by causing a 
brief electric arc to occur between the tips of 
the two conductors. Of course it would be nec- 
essary that the ends of the conductors be con- 
sumed during the reaction such that the arcs 
occur with in the slag layer. (This might be 
accomplished by using thin wires with low 
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melting points, or even using aluminized my- 
lar.) The need for an electrical power supply 
and controlling circuit is definitely a drawback 
for this method. However, the ability to have 
nearly complete control of flash timing might 
be sufficient compensation. (Note that the in- 
tensity of the light flashes would be expected 
to be roughly proportional to the period be- 
tween flashes, as this roughly corresponds to 
the amount of reactive slag produced for the 
next flash reaction.) 


Low Frequency 
Strobe Composition 


Control Signal 
Source 


Power Supply 


Figure 4. Illustration of electrically triggered 
strobe device. 


The above approaches to controlling large 
diameter strobe devices may have no merit, 
and there are considerable technical difficulties 
to overcome in designing a reliable electrical 
device. However, since the author will not be 
able to investigate this, and because the idea 
may be of use, it was thought to be appropriate 
to share the idea with other researchers. 
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Estimating the Distribution 
of Molecular Energies 


Wesley D. Smith 


Department of Chemistry, Ricks College, 
Rexburg, ID 83460-0500 


For a reactant molecule to participate in a 
chemical change, it must possess, among other 
things, the requisite activation energy E,. It 
may be useful!'! to know approximately what 
fraction of all molecules in any collection (at a 
given temperature) has at least that much en- 


ergy. 


The differential fraction of molecules that 
has a particular energy € is given by the Max- 
well—Bolzmann distribution," 


20 ef RT Je de 
(nRT)”” 


where R equals 8.314x10° kJ/mol-K, the gas 
constant, and T is the absolute temperature. So 
the fraction of molecules N that have an energy 
E, or greater is the integral of this function 
from £, to infinity. With the substitution 
x=e/RT, the integral becomes 


2 pe _ 
N=——[ e*Vx dx 
aah 
where x, = E,/RT. 


To the best of my knowledge, the integral 
(for any bounds other than zero to infinity) can 
only be done numerically, but the calculation 
is straightforward. Having done several dozen 
integrals over a range of activation energies 
and temperatures and for a number of upper 
bounds (to approximate infinity), I find that the 
results fit the equation 


E, 
RT 


InN=a-b 


Linear regression gives a = 0.9417 and b = 
0.9672. But a quick, order-of-magnitude esti- 
mation is generally all anyone needs. For that, 
a=b=1 suffices. 
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For example, suppose the activation energy 
of interest is 84 kJ/mol (about 20 kcal/mol) at 
room temperature (298 K), then 


In N= 1 — 84/(8.314x10°x298) =—33.90 


and N = 2x10'°. That is, only about a quadril- 
lionth of the molecules under these conditions 
have enough energy to react. If the temperature 
is increased to 623 K (350 °C), the fraction is 
N ~ 2x10". In other words, that rise in tem- 
perature produces one hundred million times as 
many molecules with the necessary energy. 
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Further Comment from Fred Ryan on: 
Flash Powder Output Testing: Weak 
Confinement, Issue 4, Winter 1996, 

p 5-14. [For first comments, see pp 61— 
63, Journal of Pyrotechnics, No. 5, 1997.] 


Thank you for the technical response to my 
comments on flash testing. It appears very 
clear now that the flash salutes that you were 
testing were still supersonic at 100 meters 
while the salutes that I was testing were sonic 
at 100 meters. 


This raises an interesting point. Flash sa- 
lutes produce sound waves that are initially 
supersonic but sooner or later will drop to sonic 
velocities, if not at 100 meters in your case, 
certainly at some distance probably not too 
much greater than 100 meters, since the energy 
in the wave will drop as about I/r’. The per- 
ception of “noise” by the observer can be ap- 
preciably different if the observer is located 
within the supersonic zone or outside in the 
sonic zone. Thus an interesting study presents 
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itself. Determine the zone in which supersonic 
shock waves reach the observer, and the range 
at which the shock waves have reverted to 
sonic. Get some trained observers to judge the 
quality of the sound subjectively at those two 
zones. It would be of interest to see how the 
shape of the impulse affects their judgments! 


Comment on: 

Model Rocket Engines, Theory and 
Design, Issue No. 5, Summer 1997, 
p 1-7 


I found this article extremely interesting. It 
was a pleasure to see the treatment of a cored 
endburner, which is a problem I’ve been reluc- 
tant to tackle. Now I can set up my spreadsheet 
with ease. 


I have several comments regarding the arti- 
cle. First, the figures for Example 2, Table 3 
are not totally correct. The nozzle diameter 
should be 1.70 mm and the cylindrical cavity 
height 0.37 mm, with other values changing as 
well. 


In addition to the actual values of f,,.. and 
Smin, 1 would suggest that the ratio of these two 
should be chosen with care. For “poor” propel- 
lant the ratio fira/fnin Should be kept small, 
since the burn rate exponent is large for such 
propellant. When the motor operates well at 
Jnin, there is danger of overpressurization at the 
large fnnax. The ratio can never be smaller than 
approximately 2 for an endburner, because for 
a cylindrical core of zero height, 4,,,, 1s virtu- 
ally a hemisphere of area 2nr’, while Amin 1S a 
circle of area nu’: A mad Amin=2- 


For the benefit of those using spreadsheets, 
here is an equation for the chamber pressure P, 
derived from equations 12—15 from the article: 


log N 


P.=10'” 


nd’U, «0.00001 - p,/RT. 
4bA 


noz 


where N = 


Note that neither burn time nor propellant mass 
need be known. Once the chamber pressure is 
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found, the burn rate U can be determined using 
equation 14. Finally, one fixes either propel- 
lant mass or burn time and calculates the other 
using equation 15. Constructing a spreadsheet 
for the calculations then becomes straightfor- 
ward. “Playing around” with the spreadsheet 
will provide insight into what can and cannot 
be done with endburners. 


The article has already helped in my con- 
struction of 24 mm endburners. My first three 
worked, in sharp contrast to my earlier trial- 
and-error efforts. Many thanks to Messrs. 
Clinger and Smith for their efforts. 


Sincerely, Terry McCreary, 
Associate Professor, Chemistry 


Authors’ Reply to McCreary Comment: 


Professor McCreary adds several useful in- 
sights to our article. We appreciate his com- 
ments. 


A nozzle diameter of 1.70 mm and a cylin- 
drical cavity height of 0.37 mm are correct in 
Example 2 if fini, is taken as 50. However, as 
we stated (rather obscurely) in the text, our 
results are for finin = 35. 


E. J. Klinger and Wesley D. Smith 
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Review of: 
Lecture Notes for 
Pyrotechnic Chemistry 


Paul E. Smith 


Director of Lecture Demonstrations, Purdue 
University Department of Chemistry, 1393 Brown, 
Box 177, West Lafayette, IN 47907-1393, USA 


The recent output of books providing accu- 
rate scientific discourse on pyrotechnics has 
been almost nonexistent. While there are still 
copies of many good texts on this subject, new 
editions to this extant literature have been rare. 
The publication of the two-part syllabus of the 
course “Lecture Notes for Pyrotechnic Chemis- 
try” is a welcome addition to the chemistry of 
pyrotechnics literature. (Note that earlier edi- 
tions were called: “Chemistry of Freworks” and 
“Chemistry of Pyrotechnics”.) The layout of 
this publication exists somewhere between a 
course outline and a book and has many strong 
points. 


Many great traditional pyrotechnic texts do 
not discuss or just briefly mention the topics of 
whistles and strobes. Their discussion in this 
publication is very strong and up to date. I es- 
pecially found the description of their pro- 
posed burn mechanisms helpful. The text covers 
the fundamental principles necessary for a ba- 
sic understanding of pyrotechnic chemistry. 
Readers with both a chemistry and non-chem- 
istry background will find this text a valuable 
addition to their collection of pyrotechnic lit- 
erature. 


Hazard management, in the last section of 
the text, has good charts highlighting hazard- 
ous chemical combinations. This strong section 
could only be strengthened by providing infor- 
mation about proper disposal of residue and 
waste chemicals. 


Because of the outline nature of the text, such 
things as references for some of the sources, 
page numbers, and an index did not appear in 
the earlier editions. However, edition 3.1 in- 
cludes page numbers and a nine page index. 
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Discussion of major topics was well done, 
although I would like to see more on the ef- 
fects of hygroscopicity. The presentation of 
oxidation states is very brief, and I would like 
to see this dealt with in more detail. Those 
without training in chemistry often find oxida- 
tion and reduction reactions a difficult subject. 


This volume should be of real benefit to all 
those involved in pyrotechnics and a welcome 
new addition to my pyrotechnics library. 


Review of Tom Perigrin’s 
Introductory Practical 
Pyrotechnics 


Tom Dimock 


120 Maple Ave., Ithaca, NY 14850, USA 
(e-mail: tad1@cornell.edu) 


As a pyrotechnist who has a rather public 
persona on the Internet, I am often called upon 
to recommend a book on the craft of fireworks 
to beginning pyrotechnicians. This was always 
a problem because the available literature was 
either too advanced or too dated to be really 
suitable. This gap in the literature was filled 
with the publication in 1996 of Introductory 
Practical Pyrotechnics by Tom Perigrin. 


In this small book of 200 information-filled 
pages, Mr. Perigrin more than lives up to the 
promise of the book’s name. Although he is a 
Professor of Chemistry at a leading University, 
he makes no presumptions on the readers back- 
ground—this book is readily accessible to any- 
one who can read. The book was designed to 
be the lab manual for a practical course in py- 
rotechnic chemistry and proceeds through a 
series of very specific projects beginning with 
making Black Powder by the “CIA method”, 
then progressing through a series of progres- 
sively more challenging projects up to the level 
of small aerial shells. 
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As the reader works through these projects 
they will acquire not only the knowledge of 
how to make fireworks devices, but also how 
to build much of their own tooling, and per- 
haps most important, how to incorporate ap- 
propriate safety measures into every step of 
their work. Given the inherently dangerous 
nature of pyrotechnic compositions, this strong 
emphasis on safety is more than welcome in an 
introductory text. 


If Mr. Perigrin had stopped with only the 
practical projects portion of the book, he would 
have produced an extraordinarily useful addi- 
tion to the pyrotechnic literature. Instead, he 
went on to add sections on basic pyrotechnic 
chemistry, properties of common pyrotechnic 
chemicals, a bibliography (including an excel- 
lent anti-bibliography of dangerous fringe 
texts), and several useful appendices. 


The section on pyrotechnic chemistry is 
very well presented and demonstrates Mr. Peri- 
grin’s talent for presenting complex material 
very clearly, without intimidating those with a 
weak chemistry background. References are 
provided to more rigorous texts, along with a 
strong recommendation that they be consulted 
before the reader proceeds into experimental 
formulations. 


The appendices are the part of the book that 
I keep coming back to. The first appendix “For- 


mulations” gives 49 pyrotechnic formulations 
selected from the pyrotechnic literature based 
upon their safety, efficacy, and reliability. The 
formulations are all given in a very clear pres- 
entation, along with references to their original 
publications. 


Appendices 2 and 3 give a concise reference 
to chemical nomenclature and chemical names 
and abbreviations which is very useful to those 
of us with a less than stellar chemistry educa- 
tion. Appendix 4 digests some of the most rele- 
vant parts of the BATF Orange Book [A7F— 
Explosives Law and Regulations, Bureau of 
Alcohol, Tobacco, and Firearms, Department 
of the Treasury, ATF P 5400.7 (6/90)] into a 
concise reference. Appendix 5 is a table of 
screen sizes, something that can be infuriat- 
ingly hard to find when you need it. I would 
suggest to Mr. Perigrin that a table of Black 
Powder granulations for the sporting and blast- 
ing grades would make an excellent Appen- 
dix 6. 


In summary, Mr. Perigrin has done the 
world of amateur pyrotechnics an immeasur- 
able service by providing this excellent intro- 
ductory text. It makes a great gift for the aspir- 
ing pyrotechnician and is filled with choice 
nuggets of information for the more experi- 
enced pyrotechnician. It is highly recom- 
mended. 
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